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We address the spectral properties of Mott insulators with orbital degrees of freedom, and inves- 
tigate cases where the orbital symmetry leads to Ising-like superexchange in the orbital sector. The 
paradigm of a hole propagating by its coupling to quantum fluctuations, known from the spin t-J 
model, then no longer applies. We find instead that when one of the two orbital flavors is immobile, 
as in the Falicov-Kimball model, trapped orbital polarons coexist with free hole propagation emerg- 
ing from the effective three-site hopping in the regime of large on-site Coulomb interaction U. The 
spectral functions are found analytically in this case within the retraceable path approximation in 
one and two dimensions. On the contrary, when both of the orbitals are active, as in the model for 
t2g electrons in two dimensions, we find propagating polarons with incoherent scattering dressing 
the moving hole and renormalizing the quasiparticle dispersion. Here, the spectral functions, calcu- 
lated using the self-consistent Born approximation, are anisotropic and depend on the orbital fiavor. 
Unbiased conclusions concerning the spectral properties are established by comparing the above 
results for the orbital t-J models with those obtained using the variational cluster approximation or 
exact diagonalization for the corresponding Hubbard models. The present work makes predictions 
concerning the essential features of photoemission spectra of certain fiuorides and vanadates. 



PACS numbers: 71.10.Fd, 72.10.Di, 72.80.Ga, 79.60.-i 

I. INTRODUCTION 

A hole propagates coherently in the valence band of 
a band insulator with an unrenormalized one-particle 
dispersion (i.e., determined by electronic structure cal- 
culations), whereas hole propagation in a Mott insula- 
tor is a nontrivial many-body problem. The paradigm 
here is a hole doped in the half filled one-band Hub- 
bard model that is a minimal model used to describe 
the parent compounds of high-Tc superconductors. Such 
a hole forms a defect in the antiferromagnetic (AF) back- 
ground, and its coherent propagation may appear then 
only on a strongly renormalized energy scale Naively, 
i.e., considering the Neel state induced by Ising-like spin 
interactions, one expects that a propagating hole would 
disturb the AF background and generate a string of bro- 
ken bonds, with ever increasing energy cost when the 
hole creates defects moving away from its initial posi- 
tion. This suggests hole confinement as realized already 
four decades agoi^ Nevertheless, the quantum nature of 
this problem leads to a new quality: a hole in the AF 
Mott insulator can propagate coherently on the energy 
scale J which controls AF quantum fluctuations,—'^'^ be- 
cause they heal the defects arising on the hole path. Cru- 
cial for this observation is the presence of transverse spin 
components cx {S^ SJ' + S^Sj'), responsible for quan- 
tum fluctuations in the SU(2)-symmetric Heisenberg ex- 
change interactions. 

In contrast, the orbital interactions which induce alter- 
nating orbital (AO) order are known to be more Ising- 
like^ (classical) and quantum fluctuations are then either 



substantially reducedi^ or even absent. Perhaps the most 
prominent example of robust orbital order occurs for de- 
generate Cg orbitals in the ferromagnetic (FM) planes of 
LaMnOaj^ Orbital interactions are there induced both 
by the lattice (due to cooperative Jahn- Teller effect) 
and by the superexchange^ — both of them are clas- 
sical and Ising-like, so the quantum fluctuations are to 
a large extent suppressed. Therefore, the orbital order 
is more robust than the spin order in two-dimensional 
(2D) modelSfiS and a higher hole concentration is re- 
quired to destroy it.— In the three-dimensional case off- 
diagonal (interorbital) hopping (cg orbital flavor is here 
not conservedi^) may even lead to an orbital liquid phase 
already at rather low hole dopingi^^ 

Another difference to the SU(2) spin model is that 
a hole in a ferromagnet with AO order of Cg orbitals 
can propagate coherently without introducing any string 
states. This coherent propagation arises here due to the 
Cg interorbital hopping, but is strongly renormalized by 
orbital excitations^^ However, in the t2g systems the in- 
terorbital hopping is forbidden by symmetry^i^ and the 
superexchange is purely Ising-like, so the above estab- 
lished mechanisms of coherent hole propagation in the 
regime of large on-site Coulomb repulsion U, i.e., for 
t U where t is the hopping element, are absent. One 
may then wonder whether a hole doped to a Mott insu- 
lator with the t2g AO order would then be conflnedf^^ 

The present paper is motivated by the above impor- 
tant difference between the spin physics and the orbital 
physics, especially prominent in the hole motion in the 
t2g orbital systemsJ^ In spin models, the Ising-like su- 
perexchange represents only a rather poor approxima- 
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tion for the SU(2)-symmctric Heisenberg spin exchange 
(as obtained in the t-J model by its derivation from the 
Hubbard modeU^) and could merely serve as a start- 
ing point for the full SU(2)-symmetric calculations In 
contrast, in the systems with t2g orbital degeneracy the 
Ising superexchange follows from a similar derivation as 
an accurate description of charge excitations in the sec- 
ond order of the perturbation theory in the regime of 
t ^ [/, in cases when only two orbitals are active (see 
below). Strictly speaking, in these realistic orbital sys- 
tems the Ising-like superexchange occurs in a two-band 
model when precisely one orbital flavor permits the hop- 
ping along each bond, as only then the orbital flavor can- 
not be exchanged and the pseudospin [T^TJ' + T^T^) 
operators are absent in the respective t-J model.— Such 
a situation is found not only in the above mentioned t2g 
model but also in three other distinct orbital models (see 
below). Hence, in the following paragraphs we give a 
brief overview of these four models with realistic Ising 
superexchange for which the spectral function of a sin- 
gle hole doped into the half-filled ground state will be 
studied in this paper. 

As a first example we introduce the spinless Falicov- 
Kimball (FK) model which describes itinerant d electrons 
coupled to localized / electrons. On sites occupied by an 
/ electron, the d electrons feel a strong on-site Coulomb 
repulsion U . The FK model can be solved exactly in 
infinite dimension, where it leads to a complex phase 
diagram including periodic ground states as well as a 
regime of phase separation. When the energies of two in- 
volved orbitals are degenerate (which is not the case in 4/ 
or 5/ materials), one finds that electron densities in the 
two orbitals are the same {rid = n/), and second order 
perturbation theory leads in the regime of large C/ to a 
strong-coupling model with one mobile flavor. One finds 
therefore the ground state with the AO order formed by 
sites occupied by d and / electrons on the two sublat- 
tices. The spectral density for the mobile d electrons can 
be obtained as well in onoSS, and two dimensionSj^ but 
computation of the / spectral density is quite involved 
even at infinite dimensioni^ Below we will discuss exact 
one-dimensional (ID) and approximate two-dimensional 
(2D) analytic results for the relevant FK models. 

A second and different realization of the effective low- 
energy model with Ising superexchange was proposed 
recently^ for FM planes in transition metal oxides with 
t2g orbital degeneracy: In this case, both orbital flavors 
are equivalent (the third one is inactive) and both allow 
for electronic hopping, however, each one permits the 
hopping along one axis only. Despite the ID character of 
the kinetic energy in such a model, the ground state at 
half filling has 2D AO order, stabilized by the Ising su- 
perexchange. Electron propagation, on the other hand, 
is strictly ID, so a hole replacing an electron with either 
orbital flavor may only move in one direction by the hop- 
ping t. Such a situation might be realized in Sr2V04, 
where the crystal field splits the t2g orbitals^^ and one 
finds indeed AO order— in the weakly FM planesj^ A 



different possible realization of such a model is found in 
cold-atom systems;^ with strongly anisotropic hopping 
in the spinless p-orbital Hubbard modeli^ 

The third realization concerns systems with peculiar 
Cg AO and is somewhat subtle. As described above in 
the "well-known" systems with active Cg orbitals (such 
as the perovskite manganites) the hole can propagate co- 
herently due to the interorbital hopping»^ However, one 
can identify two peculiar cases of Cg AO order where the 
interorbital hopping is strongly reduced between the or- 
bitals occupied in the ground state, and the question of 
the hole confinement in the Ising superexchange model 
is of high relevance. The first one is realized in the 
FM planes of K2CuF4f2^ or in the recently investigated 
Cs2AgF4i^ While the AO order is formed in these cases 
by Cg orbitals, crystal field stabilizes their particular lin- 
ear combinations with alternating — jy^ — Cg 
orbitals f2£ and thus suppresses the interorbital hopping 
present, for instance, in the ground state of the mangan- 
ites. In addition, the phase factors of two orbitals along 
each bond allow for the hopping only along one direc- 
tion in the plane^i^ so one arrives at a situation similar 
to that found for strongly correlated electrons in t2g or- 
bitals — it will be discussed in the Appendix „ Thus, the 
model called t2g model throughout this paper is expected 
to describe a wider class of transition metal oxides. 

Finally, another case where the AO order could in prin- 
ciple lead to the hole confinement is the situation encoun- 
tered in the ID Cg model in which the interorbital hop- 
ping and quantum fiuctuations are suppressed by sym- 
metry. Since this model is actually equivalent to the ID 
FK model discussed above, we will refer to it later sim- 
ply as to the " ID model" . Its study was stimulated by 
the experimental findings in the lightly hole-doped vana- 
dates such as Lai-a-Sr^i-VOai^ There the ground state for 
X = 0.1 is an insulating three-dimensional FM phase with 
AO order at x = 0.10. The physical origin of the Mott 
insulating state in the case of such a finite hole doping 
is not yet understood. Again, a natural question could 
be whether the peculiar stability of the Mott insulating 
state is related to the presence of the AO order. Since 
it is obvious that the hole can move in the AO state in 
the plane (due to the interorbital hopping, see above) it 
is interesting to verify whether the AO order along the 
third direction could block the hole motion. 

The four above cases are, to our knowledge, the only 
straightforward models in which the full superexchange 
is purely classical, i.e., where Ising superexchange fol- 
lows from the orbital symmetry, and is not an approx- 
imation to the Heisenberg Hamiltoniaui^ The superex- 
change models, however, have to be extended by the sec- 
ond order three-site hopping terms in each case. Such 
effective hopping terms arise in the same order of the 
perturbation theory when holes are present, and are nec- 
essary for a faithful representation of the spectral weight 
distribution in the one-particle and in the optical spec- 
tra of the underlying spin Hubbard modeli^ Therefore, 
also in the present orbital t-J models with Ising superex- 
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change similar terms are expected to play an important 
role and cannot be neglected. Moreover, while in the spin 
case such terms are in conflict with the quantum fluctu- 
ations and give thus only minor quantitative corrections 
to the coherent hole motion;^ in the orbital (pseudospin) 
models with Ising superexchange they become of crucial 
importance as they arc the only possible source for coher- 
ent carrier propagation (the Trugman loops^ with the 
hole repairing the defects on its path are here absent) 
and dictate possible coherent processesJ^ To get a bet- 
ter understanding of the balance between the coherent 
and incoherent processes in the orbital strong-coupling 
models (i.e. t-J models with three-site terms) we ana- 
lyze their spectral properties in some detail below. 

In order to arrive at a comprehensive and rather com- 
plete understanding of the elemental processes which 
accompany hole propagation in the orbital models, we 
not only combine numerical and analytical approaches 
but also we calculate the spectral properties both of the 
strong-coupling and of the respective Hubbard models. 
Hence, we first determine the Green's functions using the 
sclf-consistcnt Born approximation (SCBA) and the ana- 
lytical treatments applied to the above mentioned orbital 
strong-coupling models. This allows us to identify the 
dominant mechanism responsible for the quasiparticle 
(QP) behavior. We note that for finite doping such meth- 
ods as the slave boson approach, the path integral formal- 
ism, or the numerical approaches, are much more suit- 
able for the t- J-like models than the SCBAi^ However, 
in the one hole limit the SCBA gives reasonably good 
results which are in agreement with other methods.— 
We then compare these results to those obtained for the 
respective orbital Hubbard model. Here, in some cases 
we determine Green's functions numerically by use of ex- 
act diagonalization on small clusters, in others we use the 
variational cluster approach (VGA),— where a cluster is 
solved exactly and then embedded into a larger system. 
This variational approach is based on cluster perturba- 
tion theory developed in the last decade , '^^'^^1^° and cor- 
responds to taking the self-energy from a small cluster 
and optimizing it with respect to mean-field terms arising 
due to the AO order. The embedding via the self-energy 
approach^ allows us to include long-range (orbital) or- 
dering phenomena by optimizing a fictitious field due to 
the AO order. This method is appropriate for orbital 
Hubbard models with on-site interactions. Since the ex- 
act solution on the cluster is obtained for the full Hamil- 
tonian, it contains all potentially relevant processes like, 
e.g., the three-site hopping. 

The paper is organized as follows. The ID orbital 
model is introduced in Sec. [TTlwhereas its extension to the 
2D FK model is discussed in Sec. IIHI In both sections, 
we introduce the respective Hubbard-like Hamiltonian, 
derive from it the appropriate strong-coupling Hamilto- 
nian, and calculate analytically the hole Green's function 
for the AO state at half filling. Next, we introduce an 
exactly-solvable ID model with three-atom units along 
the chain (Sec. lIVp . called the ID 'centipede' model. 



The latter model (which was not mentioned above) serves 
merely as a didactic tool and explains the essence of 
string excitations present in the 2D model with t2g orbital 
flavors, discussed thoroughly in Sec. |Vl Here, again we 
start from the orbital Hubbard model, derive its strong- 
coupling version, and calculate the hole Green's function 
for the AO state at half filling, using two approximate 
methods described above: the SCBA and the VGA. In 
Sec. IVII we include longer-range hopping in the 2D t2g 
model (as expected in real materials) and discuss the 
main experimental implications of our study by calculat- 
ing the photoemission spectra of certain vanadates and 
fluorides. General conclusions are presented in Sec. IVIII 
The analysis is supplemented by the Appendixo where we 
derive the effective strong-coupling model for the above 
mentioned fluorides and prove that the t2g model dis- 
cussed in Sec. IVlmav indeed be applied to the hole mo- 
tion in the systems with a particular type of Cg orbital 
order. 



II. ID ORBITAL MODEL WITH ISING 
SUPEREXCHANGE 

A. Effective strong-coupling model 

As explained in Sec. HI the Ising-like superexchange 
follows if only one orbital flavor permits hopping along 
each bond, and the spins are polarized in the FM state. 
The simplest case which captures the essential features 
of the effective strong-coupling model with Ising-like su- 
perexchange follows from the ID orbital Hubbard model 

Hm = -t^{.a\°'i+i + ^■'^■) + U ^n.an^b , (2.1) 



where {b\) creates a spinless electron with orbital fla- 
vor a (6) at site i, and {7iM,7ii6} are electron density 
operators. On-site Coulomb repulsion U is the energy of 
a doubly occupied state (it arises as a linear combination 
of the Coulomb and Hund's exchange in the respective 
high-spin conflguration^) , and t is the nearest neighbor 
(NN) hopping element. Only electrons with orbital fla- 
vor a are mobile while the other ones with flavor h can- 
not hop. To simplify, we call below the a and h orbitals 
mobile and immobile ones, respectively. This situation 
corresponds to (spinless) interacting Cg electrons in the 
FM chaiu)^ or to the ID (spinless) FK model with de- 
generate orbitals. 

In the regime of large U , i.e. for i <C t/, second order 
perturbation theory leads to the effective strong-coupling 
Hamiltonian with Ising-like superexchange 

Hid ^Ht + Hj + , (2.2) 
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where 



h.c. 



(2.3) 



Ht 
Hj 



Here a tilde above a fermion operator indicates that the 
Hilbert space is restricted to unoccupied and singly oc- 
cupied sites, e.g. a\ = aj(l — nu)- The pseudospin oper- 
ators are defined as follows 



^J^fe^+i-in,n,+iV (2.4) 

i ^ ^ 



1 r 

2 ("'b ■ 



(2.6) 



and the superexchange constant J and the effective hop- 
ping parameter r are given by 



J 



4^2 
U 



(2.7) 



We introduced above the parameter r in order to dis- 
tinguish below between the terms which originate from 
the pseudospin superexchange and the hopping processes 
arising from the superexchange via the three-site terms 
that lead to the second or third neighbor effective hop- 
ping and contribute to the hole dispersion in the strong- 
coupling regime. Note that t is of the same order cx /U 
as J, so a priori these terms cannot be neglected. But 
similar as for the constrained hopping term (j2.4[) . their 
contribution is proportional to hole doping x. The ID t-J 
orbital model Titj = Ht + Hj, i.e. without the three-site 
hopping H^s, was solved exactly before^ and all excita- 
tions occurred to be dispersionless. Here we generalize 
this exact solution to the full strong-coupling Hamilto- 
nian (|2.2p including the three-site terms, and show that 
the spectral functions for both orbital flavors are then 
distinctly different. 



B. Analytic Green's functions 

We calculate below the exact Green's functions 
Ga{k,uj) and Gb{k,uj) which demonstrate whether and 
how a hole added to an a (or b) orbital may propagate 
coherently along a ID chain with the AO order. Interest- 
ingly, both functions can be determined analytically us- 
ing retraceable path approximation^ (RPA), which be- 
comes here exact because closed loops are absent in the 
ID system4^ 

An important simplification as compared with the spin 
case is the knowledge of the exact ground state |0) at 
half filling. As the Hamiltonian given by Eq. (|2.2p does 
contain then only the Ising superexchange, the Neel state 
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with a orbitals occupied on the sublattice A and b orbitals 
occupied on the sublattice B is an exact ground state. 
Here |vac) is the true vacuum state with no electrons, 
while |0) is the physical vacuum at half filling. 

We start with the Green's function for the hole doped 
in the mobile a orbital, 



lim 

(5^0 



1 



Hid - Eq + i6 



(2.9) 

where Eq is the energy of the physical vacuum at half 
filling |0), aj. is a Fourier transform of the {aj} operators 
with j G A, and the hole is created by the operator 



= v afE 



(2.10) 



with A^/2 being the number of sites in one sublattice. 
By construction, the above operator creates a hole (an- 
nihilates an electron) with momentum k on the A sub- 
lattice. After a hole is created, one finds that the state 
ak\0) in Eq. (|2.9p is an eigenstate of the Hamiltonian 
(|2.2p . The hopping (x t is blocked by the constraint in 
the Hilbert space, and the only two terms that contribute 
in this state are: (i) the superexchange term (|2.3p which 
gives the energy ^ J of two missing bonds, and (ii) the 
three-site hopping term (|2.5|) which contributes to the k 
dependence due to the processes shown in Fig. [T]^a) after 
Fourier transformation. As a result, one finds 



Gaik,u;) 



+ 2rcos(2fc) 



(2.11) 



Note that nu, = 1 in H^g , as in this case all the sites with 
j ^ B are occupied by b electrons in the ground state 
p.8p . The hole spectral fimction. 



Aa{k,u}) = ImGQ(fc,a;) 

TT 



(2.12) 



consists of a single dispersive state, shown as the mid- 
dle peak in Fig. [T](d). As expected, the hole is mobile 
thanks to the three-site terms and it propagates coher- 
ently with the unrenormalized bandwidth W = 4t. The 
result obtained here is identical with the one found us- 
ing the VGA for the corresponding Hubbard model (|2.ip 
(see also Fig. 5 of Ref. [ish . This confirms that both 
the orbital Hubbard model (|2.ip and its strong-coupling 
version with three-site terms (|2.2p arc equivalent and de- 
scribe precisely the same physics in the regime of t <C U. 

When one attempts to calculate the Green's function 
for a hole doped in the immobile 6 orbital, 



Gbik,uj) = lim ( 



''uj + HiD- Eo+iS 



one finds immediately that the state 
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FIG. 1: (Color online) Hole propagation in the ID strong- 
coupling model (|2.2|l . Two top panels show a hole doped 
into: (a) mobile a orbitals (empty boxes), and (b) immo- 
bile b orbitals (filled boxes). Solid (dashed) arrows indicate 
possible hopping processes with hopping elements t and r, 
respectively; in case of a hole added to the b orbital the latter 
process occurs only after the initial hopping by t, see panel 
(c). Panel (d) shows the exact spectral functions Aa{k, ui) and 
Ab{k,uj) of a hole added into the a orbital (middle dispersive 
feature between to = —OAt and uj = 0) and the 6 orbital (two 
side dispersionless maxima) as obtained from the ID strong- 
coupling model p.2|l . Parameters: J — OAt, t = O.lt, and 
peak broadening 5 = Q.Olt. The spectral functions obtained 
using the VGA for the ID Hubbard model ((2TJ with U = lOt 
(not shown) are identical to the exact result. 



Ht acting on iV'i^'') generates the first (normalized) state 



-ikR 



(a,_i-f a,+i)a]5^.|0), (2.15) 



with the hole delocalized to the neighboring j — 1 (j + 1) 
sites of the A sublattice, i.e., to the the left (right) from 
the initial hole position j in each Fourier component \bj) 

included in Eq. (j2.14[) . The remaining states 
with n > 2, which occur in the continued fraction expan- 
sion needed to evaluate the Green's function Gb{k, uj) (see 
below), are generated by acting (n — 2) times on ') 
with the three-site hopping term H^g- In this way one 
finds the set of symmetric states, with a superposition of 
the hole propagating forward (either to the left or to the 
right from the initial defect), i.e., along the same direc- 
tion as that given by the first hop which led to iV'fc )j 
cf. Fig. [Ijc). This structure of the basis set explains 
the absence of the k dependence in the Green's function 
for b orbitals, so we adopt the simplified notation Gb{(^) 
below. 

In the infinite basis generated by the above described 
procedure, the Hamiltonian matrix of the Hamiltonian 
([^ reads: 



[uj + J/2 V2t 

V2t UJ + 3J/4 T 

r w + J r 

r uj + J 

V 



(2.16) 



In order to calculate the relevant Green's function G(,(cl'), 
we need the (1,1) element of the inverse of this matrix. 
Due to the tridiagonal form of the Hamiltonian, this can 
be done even for an infinite Hilbert space, and we arrive 
at a continued fraction result: 



V4"' \uj + niu~Eo\ij 



2t' 




J- 



LU + J-- 



(2.17) 



where the whole self-similar part can be summed up to 
the self-energy which does not depend on fcs^ 



is not an eigenstate of the Hamiltonian Hiri- Here a 
hole is doped in each Fourier component in an occupied 
b orbital at site j in the ground state with AO order (|2.8p . 
When a hole is doped it can delocalize to its neighbors 
in the ID chain, as depicted in Fig. [TJ^b), so one has to 
introduce appropriate basis of states obtained when the 
single hole delocalizes along the ID chain. The hopping 



LU + J - 



U1 + J-- 



iU + J — S(cl') 



(2.18) 



This, together with Eq. (|2.17p . leads to a quadratic equa- 
tion for T.(lu) with two solutions: 

J:{oj) ^^{{u; + J)±^{cu + Jf - 4t2 } . (2.19) 
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The proper sign may be determined using the Green's 
function G'b(a^) obtained before*^ in the hmit of r = 0, 



2 + |J 



(2.20) 



In this Hmit the self-energy vanishes, S(ti^) = 0, and the 
Green's function has two poles at energies 



Finally, we arrive at the general result for t > 0: 



(2.21) 



Gfc(c^) 



1 



UJ - 



-J 



4^2 



(2.22) 

where the sign convention is fixed by comparing this 
result with the Green's function Gf\uj) ((2:201) : This 
implies that one has to select — (+) sign for lo < —J 
{lo > — J), respectively. 

Due to the obtained analytic structure of Gb{uj) the 
hole spectral function 



1 



ImGbiuj), 



(2.23) 



shown in Fig. [Ijd), also does not depend on k. For the 
realistic parameters with r < t it consists of two poles 
and the incoherent part centered around lu = — J. This 
latter contribution has rather low intensity and is thus 
invisible on the scale of Fig. [Ud), and the two peaks 
absorb almost the entire intensity. This result resembles 
the case of r = (|2.2ip . and might appear somewhat 
unexpected - we analyze it in the following Section. 



C. Hole confinement in a tiiree-site box 

First, we comment on the absence of the k dependence 
in the spectral function Ab{uj) (j2.23p . It sufhces to ana- 
lyze the hole in a 6 orbital at any finite value of J which 
induces the AO ground state (|2.8p . The hole can only 
move incoherently, because once it moves away from the 
initial site j [see Fig. [Ijb) and (c)], it creates a defect in 
the AO state which blocks its hopping by the three-site 
processes over the site j, see Eq. (|2.14p . Consequently, 
the hole may hop only in the other direction, i.e. away 
from the defect in the AO state, and in order to absorb 
eventually this orbital excitation it has to come back to 
its original position, retracing its path. In this way a 
forward and backward propagation along the ID chain 
interfere with each other, resulting in the fully incoher- 
ent spectrum of Fig. [T]Jd). 

Looking at the spectral function Ai,{uj) of a hole doped 
into the b orbital at finite r = O.U shown in Fig. [Ijd), 
one may be somewhat surprised that the result indicates 
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FIG. 2: (Color online) Spectral function Ab{ijj) of a hole doped 
into the b orbital in the ID model with: (a) r = 0, (b) r — 
0.5t, (c) T ^t, and (d) r = 2t. Dotted (solid) lines for J = 
(J = OAt), respectively, with broadening S = O.Olt. 



only two final states of the ID chain. These are the bond- 
ing and the antibonding state of a hole confined within 
a three-site box, and discussed in detail in Ref. ^ in 
the limit of r = 0. One finds that the two excitation en- 
ergies obtained for the present parameters, u = — 1.67t 
and uj = 1.17i, are indeed almost unchanged from those 
given by Eq. (|2.2ip at r = 0. We note that the third non- 
bonding state has a different symmetry and thus gives no 
contribution to Ab{Lu). 

Altogether, one finds that in the realistic regime of pa- 
rameters with T = J/4, the incoherent part of the spec- 
trum is extremely small and thus invisible in the scale of 
Fig. [Ijd). This implies that the hole is still practically 
trapped within the three-site box depicted on Fig. [U^b), 
in spite of the potential possibility of its delocalization by 
finite r. Only when the value of the three-site hopping 
T is considerably increased, the hole can escape from the 
three-site box, and may delocalize over the entire chain. 

A systematic evolution of the spectral function Ab{uj) 
with increasing r is depicted in Fig. O One observes 
that the incoherent spectral weight grows with increasing 
T and is already visible in between the two maxima for 
T = 0.5t. When the three-site hopping term approaches 
T = t, the spectrum changes in a qualitative way — both 
peaks are absorbed by the continuum, and the spectral 
density resembles the density of states of the ID chain 
with the NN hopping. For the extremely large effective 
hopping T ~ 2i the two peaks corresponding to the en- 
ergies given by Eq. (|2.2ip are entirely gone, and the 
spectrum corresponds to the incoherent delocalization of 
the hole over the ID chain. Note also that finite J results 
only in an overall shift of the spectra due to the energy 
cost of the hole excitation in the ordered ground state 
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III. 2D SPINLESS FALICOV-KIMBALL MODEL 



A. Effective strong-coupling model 



There are two essentially different ways to generalize 
the ID orbital Hubbard model with one passive orbital 
flavor to two dimensions in such a way that the superex- 
change remains still Ising-like. Either (i) one allows that 
the electrons with mobile flavor a can hop along all the 
bonds, i.e. in both directions in the square lattice, or 
(ii) one allows that a electrons can hop along the bonds 
parallel to the b axis, and b electrons can hop along the 
bonds parallel to the a axis. The first scenario leads to 
a special case of the 2D FK model (see below), while the 
second one describes spinless electrons in t2g orbitals of 
a FM plane and will be analyzed in Sec. |Vl 

In analogy to the ID model of Sec. |lTl the 2D FK 
model describes interacting electrons in mobile a and im- 
mobile 6 orbitals. 



where 

Ht = -t^(a|a^. +h.c.), (3.3) 

iv) 

i 

-^E(°Lb^«^Vb+i^-c-)' (3-5) 

i 
i 

-T^{a\^^nit,a,±^ + h.c)j. (3.6) 

i 

Here a and b are the unit vectors along the axes a and 6, 
while iJgP and h'^^ terms stand for linear and diagonal 
processes in the three-site effective hopping. The param- 
eters J and T are defined as in Eqs. (|2.7p . the orbital 
(pseudospin) operators are defined as in Eq. (|2.6p . 
and again the tilde above the fermion operators indicates 
that the Hilbcrt space is restricted to the unoccupied and 
singly occupied sites. 



(ij) 



a]a.f 



h.c. 



(3.1) 



Here we used the same notation as in Eq. (|2.ip . and 
(ij) are the bonds (pairs of NN sites) in the 2D lattice. 
This Hamiltonian shows complex physics^ and phase 
separation^! away from half filling. In contrast to the 
usual situation with large energy difference between / 
and d orbitals we will consider degenerate a and b or- 
bitals. Then the ground state at half filling (i.e., one 
electron per site) and large Coulomb interaction U is rel- 
atively straightforward to investigate, and one finds the 
robust AO order rather then phase separation. 



Again, we can perform second order perturbation the- 
ory in the regime of i ^ C/ as above. For the present 
square lattice there are two types of three-site terms 
— they contribute: (i) along a and b axes due to for- 
ward (linear) processes, and also (ii) connect next-nearest 
neighbor (NNN) sites along the diagonals of each plaque- 
tte in the 2D lattice, along two 90° paths. The resulting 
strong-coupling effective Hamiltonian reads 



Hfk — Ht 



H. 



(d) 
3s 1 



(3.2) 



B. Analytic Green's functions 

The Green's function Ga(k,w) for a hole in mobile 
a orbitals, defined by Eq. (|2.9p . can be calculated as 
straightforwardly as in the ID model of Sec. |TT1 and one 
finds: 



Ga{k,u;) = 



1 



+ ,/ + ek 

where the hole dispersion relation is given by 
£k ~ — 4T{(cosfc^ +cosfc,y)^ — l} . 



(3.7) 



(3.8) 



As in the ID model, see Eq. (|2.1ip . the hole propagates 
freely, resulting in a bandwidth of = 16t = 4J. In- 
deed, in the strong-coupling model (|3.2p the hopping to 
the sites occupied by b electrons is blocked by the con- 
straint. The spectral function Aa(k,i^) (|2.12p obtained 
for the hole in a orbitals consists thus of a single pole, as 
shown in Fig. [3{a). 

The dispersion relation of Eq. (|3.8p can be compared 
with the one of the lower Hubbard band obtained^i for 
the FK model (|3.ip in the regime of U ^ t, where one 
finds dispersion 




-4t(cos kx + cos ky^ 



(3.9) 



This result is the same (up to a nonsignificant constant) 
as the one obtained in the strong-coupling limit, see Eq. 
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FIG. 3: (Color online) The spectral functions for the 2D FK 
model, as obtained for mobile a orbitals (middle dispersive 
feature between lu = —1.6 and u = 0) and immobile b orbitals 
(two side dispersionless maxima): (a) with the RPA of Ref. 
for the 2D FK strong-coupling model (|3.2p with J = OAt and 
r = O.lt, and (b) by a numerical diagonalization of a 20 x 20 
cluster for the 2D FK Hubbard model ((33} with U = lOt. 
Peak broadening S = O.Olt. 



The Green's function obtained for the Hubbard- 
likc model (|3.ip is shown in Fig. ^h). While it qual- 
itatively agrees with the result derived in the strong- 
coupling limit (|3.8[) . it is here renormalized and gives 
a somewhat reduced bandwidth of the hole band. This 
indicates finite probability of double occupancies which 
hinder the three-site effective hopping processes, and re- 
duce the order parameter from its classical value as given 
in the Neel state (|2.8p . see also Sec. IV El for a similar dis- 
cussion concerning the 2D t2g orbital model. 

As in the ID model of Sec. [TTbI we use here the RPA 
to calculate the Green's function for a hole inserted into 
the immobile b orbital. However, the RPA is no longer 
exact in two dimensions, because also paths with closed 
loops are possible when electrons with one flavor are al- 
lowed to hop in both directions. We use a similar basis 
of states {iV'k"'')} f*^'' the ID calculation of Sec. Ill Bl 



to describe a single hole doped to the plane with the AO 
order. Starting from the Neel state as in Eq. (|2.8p . the 
first two states are defined as follows: 

^ = yf E ^"^"''^■^.10) ' (3-10) 



Here the first state \ip]^'') denotes the Fourier transformed 
states with hole doped into the immobile orbital at the 
initial position in the ground state of the 2D lattice, with 
AO order between two sublattices A and B, as in (|2.8p . 
see Fig. IH^a). This state may delocalize by the hopping t 
which interchanges the hole with an occupied a orbital to 
the left, right, down, or up from the initial site, resulting 
in the symmetric state with four external sites of the five- 
site polaron depicted in Fig. [DJb) — the resulting state 

(2) 

is denoted above by IV'k ) (|3.1ip . 

At this point, we have to introduce an approximation 
and we will consider only these states iV'k"'') with n > 2 
which are created on a Bethe lattice by 9 possible for- 
ward going steps from ^^). Therefore, all such states 
(with n > 2) are generated by the three-site effective hop- 
ping T, so each of them means that the hole has moved 

(2) 

forward by {n— 2) steps from the symmetric state |V'k )i 
being a linear combination of the configurations with the 
hole at one of the external sites in the five-site polaron 
(Fig. 5]). Hence, we make two approximations here, i.e., 
we assume that: (i) no closed loops occur (RPA), and 
(ii) the number of forward going steps is chosen to be 9 
which is the most probable number of possible forward- 
going three-site steps on a square lattice with the AO 
order.— Let us emphasize, however, that the closed loops 
which are neglected here do not lead to the delocaliza- 
tion of the hole as these are not the so-called Trugman 
loopa^^ where the hole could repair the defects by a circu- 
lar motion around a plaquette in the square lattice. The 
Hamiltonian [Eq. (j3.2p ] matrix written in this basis is as 
follows: 



4 


- Hfk 




) = 




LJ + 


J 


2t 








2t 






3r 










3r 


LJ + 2J 


3t .. 










3t 


UJ + 2J .. 



(3.12) 



Again, as in Eq. (|2.16p . the k-dependence is absent, and 
the Green's function Gb{(^) can be calculated using con- 
tinued fraction method in a similar way as in the ID case 
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FIG. 4: (Color online) Schematic view of a five-site polaron 
which occurs after the removal of one fe-electron in the 2D 
FK model (|3.2|l . (a) Initial state with a hole (empty circle) 
on B sublattice — crosses represent immobile 6 electrons and 
filled circles mobile a electrons. Five-site polaron is indicated 
by the dotted line, (b) The hole may delocalize within it to 
the nearest neighbors by hopping t (solid line), generating 
three broken bonds which cost energy 3 J/4 for each hopping 
process. The hole can hop directly between the external sites 
(dashed line) by second order three-site hopping r = J/4 via 
a doubly occupied site (outside the polaron). 



(cf. Sec. Ill B|) and we obtain 

Gb{L0) = {iJ + J 



8f 



UJ - 




where we select — (+) sign for ui < — 2J {cu > —2 J), 
respectively. 

The spectral function of a hole doped into the b orbital 
is shown in Fig. ^a). It consists of two distinct disper- 
sionless peaks, and a dispersionless incoherent part with 
negligible spectral weight [invisible on the scale of Fig. 
[3lja)] between them. While the incoherent spectrum dis- 
appears in the limit of vanishing three-site hopping r, 
the two dispersionless peaks survive and the distance be- 
tween them becomes At for J — *■ 0. Hence, the doped 
hole is not only immobile but also trapped within the 
five-site orbital polaron depicted in Fig. 4(a) - only four 



sites can be reached from the central site by NN hopping 
t, so the ground state can be found in the truncated basis 
{|V'k^^)j li'k^}} (see below). This situation resembles very 
much the ID case discussed previously, and indeed the 
same discussion as the one for the hole confinement in a 
three-site box presented in Sec. Ill CI applies here. In the 
next Section we discuss in detail the quantitative argu- 
ments which suggest the hole confinement in this five-site 
polaron. 

We would like to emphasize that the spectral functions 
obtained using the RPA for the strong-coupling model 
Eq. (|3.2p are almost identical to the ones obtained nu- 
merically by exact diagonalization of the Hubbard model 
dnH) on a 20 X 20 lattice, cf. Fig. ^h). This means 
that the crude assumption of walks without closed loops 
made within the RPA approximation (i.e., replacement 
of the square lattice by the Bethe lattice) is a posteriori 



well justified for the strong-coupling model defined by 
Eq. (|3.2p . We provide also more arguments which com- 
plete our understanding of this result in the next Section. 
Furthermore, this means that not only the RPA method 
is correct but also the two models (Hubbard and the 
strong-coupling model) are fully equivalent and describe 
the same physics in the considered regime of parameters. 

Finally, we note that the results obtained by the VGA 
(not shown) are very similar to the ones of exact diago- 
nalization, if we use periodic boundary conditions. While 
open boundaries are usually optimal for the VCA^ they 
can 'cut' the five-site polaron and lead to signals at 
wrong frequencies. In a large enough cluster, these con- 
tributions from polarons with less than five sites would 
have vanishing weight, but for the cluster sizes consid- 
ered here, self-energies with periodic boundary conditions 
have to be used to eliminate them. 



C. Localized five-site orbital polaron 

The following comparison shows that the two domi- 
nant peaks of the Green's function for b orbitals Gb{uj) 
can be well reproduced by taking into account the po- 
laron alone, i.e., by considering just a cluster of five sites 
depicted in Fig. |4(b)] and neglecting the rest of the lattice. 
In this case the problem can be solved by diagonalizing 
the Hamiltonian in the basis consisting of two states de- 
fined in the last section: |?/'^^^) (PUT)) and (PTT|) . 
This means that the infinite matrix of Eq. (|3.12p for 
the hole doped into the central site of the polaron re- 
duces to the 2x2 matrix and one obtains the energies of 
two poles of the Green's function Gb{i-^), corresponding 
to the bonding and antibonding state within the five-site 
polaron: 
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^1,2 



-J-T±\ Ut 



9 

64' 



J2 



(3.14) 



Assuming J = OAt and r = O.lt in Eq. (|3.14p we 
obtain uji_2 = {1.37t, — 2.67t}. This compares very favor- 
ably with the results obtained for the strong-coupling 
model (I3I2D: (i) within the RPA (see Sec. UTTB]) . 
tJi 2 = {1.38t, — 2.69i}, and (ii) using the numerical anal- 
ysis of this model on the 20 x 20 lattice, which gives 
^1,2 = {1.40^, -2. 74t} (not shown). We stress that the 
excellent agreement between all these methods demon- 
strates that the RPA (i.e., full continued fraction) turns 
out to be only slightly better than the calculation re- 
stricted to the five-site polaron of Fig. 4(b) It means 



that the probability of the configurations with the hole 
outside the five-site polaron is indeed very low, and it 
explains why the RPA assumption of having no walks 
with closed loops works here so well. Lastly, we note 
that all these results agree quite well with the numerical 
ones for the itinerant FK model (|3.ip . cf. Fig. [3jb) with 
the peaks situated at lox = lA2t and L02 = — 2.56t. 

The present five-site orbital polaron resembles the five- 
site spin polaron identified in Monte Garlo studies for the 



10 



2D Kondo model4^ For example, as for the spin polaron 
in the Kondo model, the spectral density of the orbital 
polaron is comprised of two dispersionless peaks with a 
distance of At for J ^ 0. There is, however, one differ- 
ence: Here not only the hole can move by direct NN hop- 
ping t between the central site and the four external sites, 
but there is also a second order three-site diagonal hop- 
ping (j3.6p which couples directly the neighboring external 
sites of the polaron (see Fig.|4(b)|). and contributes to the 

(2) 

energy of the \tpl^ ) state. Actually, due to the inclusion 
of these processes (which enable the smallest loops on 
the lattice, with two t and one r hopping processes) we 
could obtain the above mentioned perfect agreement be- 
tween the numerical, the RPA, and the five-site polaron 
results for the strong-coupling version of the FK model 
Eq. p.2p . Otherwise, e.g. for J = OAt the energies of 
the two peaks in the RPA (five-site polaron) calculation 
would be equal to {1.47^,-2.580 [{1.46^,-2.56*}], re- 
spectively, and would only rather poorly agree with the 
numerical results of Eq. (|3.2p . 

Summarizing, the holes doped into the immobile or- 
bitals of the FK model are almost entirely localized 
within the five-site orbital polaron of Fig. 4(b) In order 



to calculate the energy of this polaron correctly one has 
to take into account the energies of the processes which 
involve four external polaron sites. In addition, we note 
that the widely used SCBA^ (used for the orbital strong- 
coupling model in the next Section) does not work so 
well for the case of hole doped into the immobile orbital 
of the FK model, as it does not respect the constraint on 
the hole motion. It incorrectly uses an on-site energy of 
J instead of 3 J/4 for the excitations at external sites of 
the polaron. As shown above, the hole spends almost all 
its time inside the polaron and hence this underestima- 
tion of the energy heavily influences the energies of the 
poles of the Green's function in this case [e.g. the lowest 
peak for J = OAt is situated almost at —3t in the SCBA 
calculations (not shown)]. 



IV. STRING EXCITATIONS IN THE ID MODEL 



to generate some characteristic features recognized later 
in the spectral properties of the 2D t2g model (see Sec. 
EJ. 

The ID model of Fig. [Sja) consists of a chain along b 
axis, with the Hamiltonian as described by Eq. (j2.ip . and 
two sites being the NNs of every second site of the chain 
along the a axis, which could represent radicals added to 
a linear molecule. We use here the convention introduced 
before for the t2g orbital systems^^i^ that a and b or- 
bitals stand for yz and zx t2g orbitals, respectively, that 
permit the electron hopping along the b and a axis in the 
(a, b) plane. The Hamiltonian of the present (called here 
centipede) model is, 

= -*E{4(^2vu+^2..d)+h.C.} 

i 

- t^^{a\a^^i+\\.c.) -VU^^niaUib ■ (4.1) 

i i 

The hopping along the bonds parallel to the a axis is 
allowed only to the orbitals &, with the corresponding 

creation operators {&2i „, &2i d}- ^'^'^ ^^S- 151 s-)- To simplify 
notation, we call these orbitals u and d, and introduce the 
following operators: 



"2i.d ■ 



(4.2) 



In the limit of large U {U ^ t) the occupied orbitals 
form AO order along the chain and we select the Neel 
state with b {u and d) orbitals occupied on the external 
sites, as shown in Fig. [5l as we are interested in their 
effect on the hole propagation when it was doped to an 
a orbital. This leads to the following strong-coupling 
version of the ID centipede model (|4.ip : 

i 

- T5I("2i^2j+l,&a2,+2 + h.c.) 



(4.3) 



The ID model (Sec. H]) and the 2D FK model (Sec. 
imp bear the same generic features: (i) a hole generated 
in the so-called mobile orbital always leads to the disper- 
sive spectrum with the full unrenormalized bandwidth, 
and (ii) a hole doped into the so-called immobile orbital 
is localized, leading to a non-dispersive spectral function. 
In this Section we investigate the consequences of string 
excitations which may arise when both orbital states al- 
low only ID hopping, as in the case of two t2g orbitals 
lying in two vertical planes with respect to the considered 
plane. Thus we will study the ID model with electrons 
hopping between yz and zx orbitals in (a, b) plane — 
the model has only 2N sites for the chain of length N, 
see Fig. [5j We will show that even the shortest possible 
strings with the length of one bond which can be excited 
here when the hole moves in this geometry are sufficient 



On the one hand, the superexchange interaction for all 
the bonds within the centipede was not included in Eq. 
(|4.3p as it results only in a rather trivial energy shift of 
the spectra obtained from the Green's function Ga{k,uj) 
which is of interest here^^ cf. Sec. Ill Bl On the other 
hand, the last term in Eq. (|4.3p was added to simulate 
the creation of string excitations which occur in the full 
2D model of Sec. |V] (see also discussion below). 

Whereas the second term in Eq. (|4.3p is once again 
the three-site hopping derived before in the ID model 
(|2.5p [cf. Fig. [5ja)], the other two terms describe the 
possibility of creating defects in the AO order when the 
hole leaves the spine of the centipede (i.e., moves away 
from the a sites) by creating strings of length one, just as 
it may happen in the t2g 2D model, sec Sec. |Vl Here the 
hole can leave the chain to its NN orbital U2i or d2i [cf. 
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FIG. 5: (Color online) Propagation of a hole added into the 
a orbital in the centipede strong-coupling model (|4.3p : (a) 
schematic picture of a hole doped at site a and its possible 
derealization via hopping t (solid lines) and three-site ef- 
fective T term (dashed lines); (b) spectral function Aa{k,u!). 
Parameters: J = OAt, r = O.lt, peak broadening S = O.Olt. 
The chain is oriented along the b axis, and nonequivalent po- 
sitions of the orbitals which do not permit hopping along this 
direction are labelled b, u and d in panel (a). 



sites attached to the chain along the a axis shown in Fig. 
[5]Ja)]. Such defects are created by hopping t and costs 
energy 3J/4 in each case. Hence, the present ID model 
represents an extreme reduction of the full <2g 2D model, 
allowing only the strings of length one, and each defect 
has to be deexcited before the hole can hop to another 
three-site unit along the chain. Note however, that the 
energies of these string excitations are properly chosen 
and are just the same as in the full 2D model of Sec. |Vl 
The model given by Eq. (|4.3p constitutes a one-particle 
problem (after inserting n2i+i^b = 1 which is consistent 
with the Ising nature of the superexchange) and hence 
can be solved exactly. We will consider the Green's func- 
tion Ga{k,Lo) for a orbitals, defined similarly as in Eq. 
(|2.9p , and a hole excitation is created again by the oper- 
ator ak of Eq. (|2.10p . The continued fraction terminates 
after the second step and one finds the exact Green's 
function 



Ga{k,i0) 



1 



2TC0s(2fc) 



2*2 



(4.4) 



leading to the corresponding spectral function ^^(fc,^), 
defined as in Eq. (|2.12p . The numerical results obtained 
with J = 0.4t are shown in Fig. [IJb). Instead of a single 




FIG. 6: (Color online) Characteristic features in the spectra 
obtained for the ID centipede model (Fig. [5]) for increasing 
J: (a) the bandwidth Wi,2, (b) the spectral weight aqp, and 
(c) the distance A between the two peaks. The solid (dotted) 
line in (a) corresponds to the first (second) dispersive peak in 
Aa{k,uj) whereas the solid (dashed) lines in the lower panels 
show results for fc = (fc = 7i"/2), respectively. The light 
solid line in (a) is merely a guideline for the eye to show that 
the bandwidth of the first peak is a function with a positive 
second derivative. Parameter: t = J/4. 



dispersive state of Fig. [TJd), the spectral function con- 
sists here of two dispersive peaks, separated by a gap of 
roughly 2\/2t. This demonstrates that the larger hopping 
t suppresses at first instance the hopping along the chain 
by the element t, and a hole doped into the a orbital de- 
localizes in first place over the three-site unit, discussed 
in Sec. Ill CI consisting of a hole and two b {u and d) or- 
bitals. Therefore, the hole behaves effectively as a defect 
created at a 6 site in the ID chain of Sec. [Ill This ex- 
plains that the maxima of Aa(fc, uj) are found again for a 
bonding and antibonding state, similar to the structure of 
Ab{uj) in Sec. Ill Bl However, at present the correspond- 
ing states gain weak dispersion because the hole may as 
well delocalize along the chain by the three-site hopping 
r. Note also that the low-energy (right) peak has slightly 
higher dispersion (leading to a broader band) than the 
left one. This case illustrates that the ID dispersion is 
broader for the QP state but is also shared by the fea- 
ture at higher energy. This observation will help us to 
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interpret the spectra for the 2D t2g model in Sec. El 

In addition we also calculated some characteristic fea- 
tures of the spectra of the centipede model, cf. Fig. [6l 
They will mostly serve for a comparison with the respec- 
tive results of the 2D t2g model, presented in Sec. IV El 
However, let us only remark that the renormalization of 
the bandwidth, shown in Fig. [Hl^a) follows from an in- 
tricate interplay between coherent hole propagation and 
the string excitations. With increasing r = J/4 the free 
bandwidth increases but at the same time the energies 
of the defects (generated by the hole when it moves to 
'lower' or 'upper' sites) are cx J; hence, the bandwidth 
does not depend in a linear way on J, cf. Fig. ^a). 
Physically this means that the hole motion is gradually 
more and more confined to just the ID path along the 
chain with increasing J (and keeping r = J/4) . 

V. 2D MODEL FOR t2g ELECTRONS 
A. Effective strong-coupling model 

After analyzing the spectral properties of the simpler 
ID model and 2D FK model, we consider below the model 
relevant for transition metal oxides with active t2g or- 
bitals, when the crystal field splits them into Cg and oi 
states, and the doublet Cg is filled by one electron per 
site. This occurs for the configuration (e.g. in the 
titanates) when the eg doublet has lower energy than 
the ai state, or for (P configuration when the Cg states 
have higher energy and are considered here, while the ai 
state is occupied by one electron at each site and thus 
inactive (as in the high-spin ground state of the i?V03 
perovskiteS)^ where R stands for a rare earth element). 
To be specific, we consider electrons with two t2g orbital 
flavors, yz = a and zx = b, moving within the (a, b) 
plane. In contrast to the FK model with two nonequiv- 
alent orbitals and only one orbital flavor contributing to 
the kinetic energy (Sec. IIIip . both t2g orbitals are here 
equivalent and electrons can propagate conserving the 
orbital flavor by the NN hopping t, but only along one 
direction in the (a, b) plane4^ While this results is a com- 
plicated many-body problem at arbitrary electron flUing, 
the motion of a single hole added at half filling remains 
stiU strictly IDJ^ 

The orbital Hubbard model for spinless electrons in the 
FM {a,b) plane reads: 

(■y>l|a (U>l|b 

+ U^niaUib, (5.1) 

i 

where a and b are the orbital flavors with the same hop- 
ping t along b and a axis, respectively, and U stands again 
for the on-site interaction energy for a doubly occupied 
configuration. At the filling of one electron in {a, 6} or- 
bitals per site this interaction corresponds to the high- 



spin (P (or (P) state. Second order perturbation theory 
applied to this Hamiltonian in the regime oit leads 
then to the strong-coupling model, 

nt,^^Ht+Hj + H^,'J +Hi^\ (5.2) 

where 

Ht = -iE(^I^«+a + «I«.+b+^^-C-), (5.3) 

i 

i 

~'^E("Lb"'*''"i+b + ^^■'^■) ' (5-^) 

i 

= -^E(«Lb^''^^^^±a + h-C.) 

i 

-^E(«U^^^^^^±a + l^-^-)- (5.6) 

i 

The parameters J and r are defined as in Eqs. (|2.7p . 
whereas the pseudospin operators Tf are defined as in 
Eq. (|2.6p . Again the tilde above the fermion operators 
indicates that the Hilbert space is restricted to the unoc- 
cupied and singly occupied sites. One interesting obser- 
vation here is that the strictly ID kinetic energy of the 
two orbitals leads to the 2D superexchange (|5.4p . As in 
the spin case?^ the superexchange is active only when 
electrons with two different flavors occupy the neighbor- 
ing sites (one bond), but here only one of them can hop, 
which explains the prefactor i in Eq. (|5.4p . 

Instead of the quantum behavior and frustration 
present in the compass modelfS here one flnds that the 
perfect AO ordered state |0) (|2.8p is the ground state of 
the model at half flUing. Figure [7] presents in a schematic 
way a few flrst steps in the motion of a hole inserted at 
a selected site into such a ground state. When the hole 
moves via NN hopping t, it creates string excitations in 
each step that cannot be healed by orbital flips, because 
the orbital superexchange (|5.4p is purely Ising-like, see 
Fig. Efc). Moreover, it cannot heal the defects by itself, 
because it cannot complete a Trugman loop^ when the 
orbital defects are created and three occupied orbitals are 
moved anticlockwise on a plaquette after the hole moved 
clockwise by three steps, see Fig. [DJd). 

The structure of Titj^ (|5.2p is similar to that of the 2D 
FK model — one obtains again three-site terms in the 
strong-coupling model along the axes (to third neighbors) 
()5.5p as well as along the plaquette diagonals (to second 
neighbors) (j5.6p . However, there is now an important dif- 
ference to the 2D FK model: While three-site terms along 
the axes, given in Eq. (j5.5p . involve a single electron and 
conserve orbital flavor, terms along the diagonal (j5.6p 
require the subsequent hopping of two electrons with dif- 
ferent orbital flavor in each step, so the orbital flavor (at 
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FIG. 7: (Color online) Schematic view of the hole motion 
in the strong coupling t2g orbital model (jS-lf) with AO or- 
der. Circles depict holes while horizontal (vertical) rectangles 
depict occupied orbitals with electrons that can move only 
horizontally (vertically), respectively. The hole inserted in 
the AO state (a) can move via NN hopping t, but it has to 
turn by 90° in each step along its path and leaves behind bro- 
ken bonds leading to string excitations with ever increasing 
energy (b) and (c). After moving by 270° around a plaquette 
(d), the hole cannot return to its initial position as would be 
necessary to complete the Trugman path.— 



the site where the double occupancy is created in the ex- 
cited state) is flipped. We will see later that such terms 
flipping orbital flavor arc in fact suppressed, because, in 
contrast to the forward hopping along one of the cubic 
axes, they disturb the AO order in the background and 
thus cost energy. As such, they do not affect the low- 
energy QP state, but contribute only to the incoherent 
processes at higher energy. 

To achieve a complete understanding of the excitation 
spectra at half filling, we used two complementary meth- 
ods and investigated the orbital Hubbard model (|5.ip us- 
ing the VGA, and the strong-coupling model (|5.2p within 
the SCBA. Both cases were supplemented by exact diag- 
onalization for small (4x4 and 4x6 sites) clusters. In 
the following Section we formulate the SCBA treatment 
of the t2g model, next give results of the SCBA calcu- 
lations in Sec. IV C[ compare them to a numerical VCA 
treatment in Sec. IV Dl and discuss the QP properties in 
Sec. IV El Finally, the impact of longer-range hopping 
pertinent to realistic materials is treated in Sec. IVII 



B. Hole-orbiton coupling in the t2g model 

The calculation of spectral properties of the strong- 
coupling model given by Eq. (|5.2p is more involved than 
in the previous cases. On the one hand, in each such step 
by hopping t the position of a hole is interchanged with an 
electron and a defect in the AO state is created (see Fig. 
[7]). Therefore, one arrives at a situation analogous to a 



hole which tries to propagate in an antiferromagnet with 
the Ising interactionsi^ On the other hand, the important 
new feature which makes the present t2g problem more 
complex is that the electron hopping t is now allowed 
for both orbital flavors. Thus, this problem cannot be 
solved by the RPA^, which was used to determine the 
Green's functions Ga(ii.,uj) and G{,(k, w) of a hole doped 
into the a and b orbital in cither the ID models or in 
the 2D FK model. Moreover, this problem cannot be 
reduced to any effective one-particle Hamiltonian that 
one could solve at least numerically for large clusters. 
Hence, we use below the SCBA which gives quite reliable 
results in the spin casci^ Here one finds that it is well 
designed to treat this problem because several processes 
not included in the SCBA drop out of the Hamiltonian 
(|5.2p for physical reasons (sec below), and therefore the 
approximation performs remarkably well. 

In order to implement the SCBA we have to reduce 
the model of Eq. (|5.2p into the polaron problem, follow- 
ing Ref. H. Firstly, we divide the square lattice into two 
sublattices A and B, such that all the a {b) orbitals are 
occupied in the perfect AO state in sublattice A {B), re- 
spectively, see Eq. (|2.8p . Secondly, we rotate the orbital 
pscudospins on the A sublattice [corresponding the down 
orbital fiavor, see Eq. (|2.6p ] so that all the pseudospin 
operators take a positive value, (T/) = 1/2, in the trans- 
formed ground state. Finally, we introduce boson oper- 
ators ai (responsible for orbital excitations - orbitona^) 
and fermion operators hi (holons), which are related to 
the ones in the original Hilbcrt space by the following 
transformation: 



hlil-aja,), 



a. = h]a,. 



(5.7) 



Note, that we added the projection operators (1 — aja.j) 
to the transformation relation for the bi fermions in order 
to keep track of the violation of the local constraint that 
'no hole and orbiton can be present at the same site', cf. 
constraint CI in Ref. [^. 

Before writing down the polaronic Hamiltonian, we 
make the following approximations: (i) keep only linear 
terms in boson operators (as we use linear orbital-wave 
approximation^), (ii) skip (1 — hi hi) and {1 — alai) pro- 
jection operators when deriving the effective Hamiltonian 
(both simplifications are allowed for the present case of 
one hole and Ising superexchangc) , and (iii) neglect the 
orbital-flipping terms Eq. (j5.6p as generating the cou- 
pling between the hole and two orbitons and leading to 
higher-order processes in the perturbation theory. Then, 
after Fourier transformation, the Hamiltonian Eq. (|5.2p 
reads. 



He 



Ht + Hj + H. 



(0 

3s ' 



(5.8) 
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with 



Ht 



z 
N 



X] |^^(k,q)/ik^/ik-qBaqA 

k,q 

-iV(k, q)4B/ik-qAaqB + li.c.} , (5.9) 



Hj = ujQ 



E 



(5.10) 



^3s = E{-^-4(k)/iL/lkA + eB(k)4s/ikB},(5.11) 



where z = 4 is the coordination number of the square 
lattice, the sums are over momenta k in the fuU Brillouin 
zone for the whole lattice^^ the total number of sites in 
the plane is A'', and indices A and B denote the orbiton 
operators in both sublatticcs. The orbiton energy loq = J 
docs not depend on momentum k, and the vertices in Eq. 
()5.9p have ID dependence on momenta: 



M(k,q) = itcos(fc^ - q^) , 

-^(k,q) = ^tCOS(fc.y - Qy) , 



(5.12) 
(5.13) 



whereas the ID hole dispersion arising from the propa- 
gation within the sublatticcs in Eq. (|5.1ip are: 

£A(k) = 2TCOs(2fcy) , (5.14) 
£i3(k) = 2TCOs(2fc^) . (5.15) 

C. Self-consistent Born approximation 

Instead of calculating hole Green's functions Ga(k, w) 
and Gb(k, cj) using their definitions (see Sec. |TT|, it is 
convenient now to express them in terms of the operators 
used in Eq. (|5.8p . Hence, we introduce hole creation 
operators on sublattice B 



h 



kS 



-ikRi i,t 



(5.16) 



Next, using Eqs. (|5.7p we obtain the relation: 



«=/|E' 

jeB 



ikR 



=yiE^""'''^^'(i-"i".)io)=4Bio), 

(5.17) 

since one does not have any pseudospin defects in the AO 
ordered state |0) (|2.8p . The latter feature is also respon- 
sible for the fact that one cannot annihilate an electron 
with the 'wrong' flavor, e.g. in the b orbital on the A sub- 
lattice in the ground state |0), which justifies the above 
definition of the Fourier transformation. While one still 
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FIG. 8; Diagrammatic representation of the perturbative pro- 
cedure used within the SCBA: top — the Dyson's equation 
for the Gbb{^,^^) and Gaa (k,ti^) Green's functions; bottom 
— the summation of diagrams for the self-energy E_B_B(k, w). 
The densely-dotted and the dashed-dotted rainbow lines in 
the self-energy (lower part) connect the two vertices A''(k, q) 
and Af(k, q), respectively. 



needs to perform rotation of the pseudospin flavor on 
sublattice A, a similar relation can be obtained for 
operators. Finally, we obtain that 



GAA(k,c^) = lim(0 

(5-^0 



GBB(k,cj) = lim(0 



^kA 



'kS 



fat 

1 





(5.18) 



■h 





(5.19) 



We calculate the above Green's functions (|5.18p and 
(|5.19p by summing over all possible noncrossing diagrams 
(i.e., neglecting closed loops), cf. lower part of Fig. [H 
However, the crossing diagrams do not contribute here 
since the closed loops (Trugman processes) do not occur, 
see Fig. [71 Since the structure of the present problem 
makes it necessary that two Green's functions and two 
self-energies are considered, we write the Dyson's equa- 
tion for each of them, as represented in Fig. [8] 

G^^(k,c.) = {G^"l(k,c.)}"'-I]AA(k,c.),(5.20) 
GB^(k,c^) = {G<5"]5(k,c^)}"'-EBB(k,c^),(5.21) 

where the free Green's functions are given by. 



G^°i(k,c.) 



1 



+ J + e^(k) ' 
1 

w + J + es(k) ' 



(5.22) 
(5.23) 
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and the self-energies 

7? 

EAA(k,w) = — ^M2(k,q) Gi3B(k-q,w-c^o), 
q 

(5.24) 

SBB(k,c^) = ^^iV2(k,q) G^A(k-q,c^-c^o)- 
q 

(5.25) 

are obtained by summing up the rainbow diagrams of 
Fig. [HI Note that the intersublattice Green's function 
G AB (k, io) vanishes since it would imply that at least 
one defect was left in the sublattice B after the hole was 
annihilated in the sublattice A, resulting in orthogonal 
states as there arc no processes in the Hamiltonian which 
cure such defects [cf. the form of the Hamiltonian Eq. 
((5:81) and Fig. [8]. 

We solved Eqs. ([ODl) - ([OT|) together with Eqs. 
(|5?24l) - ([5:25l) self-consistently on a mesh of 20 x 20 k- 
points (and checked the convergence comparing the re- 
sults with those obtained for the cluster with 32 x 32 
k-points). The spectral functions defined for the sublat- 
tices 

A,(k,c^) = -ilmG^^(k,u;), (5.26) 
A,(k,c^) = -ilmGsi3(k,a;), (5.27) 

TT 

are displayed in Fig. (5] As discussed in detail in Ref. 
[isl . the spectral density consists of dispersive ladder-like 
spectrum suggesting that the hole doped into any of the 
two orbitals is mobile. The dispersion is particularly pro- 
nounced for the first (low-energy) excitation which we 
identify as a QP state. One finds that its dispersion is 
strictly ID and is dictated by the orbital fiavor at the site 
where the hole was added, i.e. no dispersion occurs in the 
complementary direction. For example, a hole added to 
the a orbital moves (thanks to the three-site terms) only 
along the h direction. However, such a hole moving along 
the h direction due to the three-site terms could also un- 
dergo incoherent scattering on orbital excitations, and in 
addition performs "excursions" to the B sublattice due 
to the t processes, which create string-like states. The 
peculiar interrelation of these two types of (coherent and 
incoherent) propagation (which we discuss in detail in 
Sec. IV Ep leads to the spectra depicted in Fig. [Q] 

The lack of the QP dispersion in one direction, e.g. 
along the a direction for a hole doped into the a orbital, is 
at first instance counterintuitive: One could imagine that 
it should be allowed that the hole doped into the a orbital 
switches to a neighboring site of the B sublattice by the 
t process, and then propagates freely along the a axis by 
the three-site effective hopping t without generating any 
further defects. This might lead to some dispersion in the 
spectra along the k^, direction. However, the hole always 
has to return to the original site where it has been doped 
as it has to erase the defect it created in the first t step 
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FIG. 9: Spectral function as obtained in the SCBA for the 
effective t2g model (|5.8p for a hole doped into: (a) a orbital, 
and (b) h orbital. Parameters: J = 0.4t, r = O.lt, and peak 
broadening 5 = 0.0 If. 

when it moved to the other sublattice (otherwise, the hole 
annihilation operator would not permit to return to the 
ground state) . Note that this behaviour is similar to the 
hole confinement in a three-site cluster as calculated for 
the hole doped into the 6 orbital in the ID model, cf. Fig. 
[T}^ As a result of such processes, one finds very small 
incoherent (and k- independent) spectral weight in the 
spectra of Fig. [9l which remains invisible in the present 
scale. This discussion demonstrates also that the spectra 
found for the 2D tig orbital model are dominated by the 
ID physics explained in Sees. HIl and [TVl 

Next, three remarks which concern the validity of our 
results are in order here. Firstly, note that if we skip the 
flavor-conserving three-site terms (|5.11l) . the calculated 
spectral functions (not shown) reproduce the well-known 
ladder spectra and arc equivalent to those calculated for 
the Ising limit of the spin i-J modeli^ This means that 
the zig-zag-like hole trapping in the orbital case is phys- 
ically similar to the standard hole trapping in the spin 
case (apart from the modified energy scale due to a dif- 
ferent value of the superexchange, the ladder spectra are 
similar in both cases), whereas for the free hole move- 
ment obviously it matters whether the dispersion rela- 
tion is ID or 2D. Moreover, this also means that in this 
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FIG. 10: (Color online) Schematic representation of two 
three-site terms in the t2g orbital model (|5.2|) . Circles de- 
pict holes while horizontal (vertical) rectangles depict occu- 
pied orbitals with electrons that can move only horizontally 
(vertically), respectively. Processes shown in panels (a)-(c) 
result from forward propagation (|5.5|l . while the ones shown 
in panels (d)-(f) and given by Eq.. (|5.6|) create a defect in the 
AO order with the energy cost indicated by the lines (broken 
bonds) in (f). 



special case (r = 0) the spectra are the same for holes 
doped into either of the orbitals as the Green's functions 
are the same for both sublattices. However, even in this 
case it is not allowed to assume a priori that A = B 
and Gy!iA(k, cj) = Gbb{^iI^)- In fact, these are two sub- 
lattices with two distinct orbital states occupied in the 
ground state at half filling, and each orbital has entirely 
different hopping geometry. This does not happen in the 
standard spin case with isotropic hopping, and for this 
reason one can eliminate there the sublatticc indices. 

Secondly, to obtain the result shown in Fig. [5] we ne- 
glected the three-site terms with 90° hopping, see Eq. 
(|5.6|) . One may wonder whether this approximation is 
justified whereas the formally quite similar forward hop- 
ping term (|5.5p is crucial and is responsible for the ab- 
sence of hole confinement in the ground state with the 
AO orderJ^ Hence, let us look in more detail at these 
two different kinds of three-site terms, shown in Fig. [TOl 
The first (linear) hopping term (|5.5p transports an a elec- 
tron along the b axis over a site occupied by a 6 electron. 
Such processes are responsible for the ID coherent hole 
propagation. As we can see in Figs. fTUT a^-fc). the AO 
order remains then undisturbed, so these processes deter- 
mine the low-energy features in the spectra. Hopping by 
the other three-site term (|5.6p . shown in Fig. [TOld)-(f). 
involves an orbital fiip at the intermediate site, destroys 
the AO order on six neighboring bonds, and thus costs 
additional energy. As two orbitals are fiipped and two ex- 



cited states are generated, these processes go beyond the 
lowest order perturbation theory, and it is consistent to 
neglect them in the SCBA. In any case, they could con- 
tribute only to the incoherent processes at high energy 
and not to the low-energy QP. Indeed, this interpretation 
was confirmed by exact diagonalization performed for the 
strong-coupling Hamiltonian (|5.2p on 4 x 4 and 4x6 clus- 
ters, which gave the same results for the QP dispersion, 
no matter whether the orbital-flipping terms (|5.6p were 
included or not. In addition, the QP dispersion found in 
the SCBA agrees with the numerical results obtained by 
the VGA (see below), which gives further support to the 
present SGBA. 

Lastly, despite several other approximations made in 
writing down the Hamiltonian Eq. (|5.8p . the vertex part 
Ht is exact, in contrast to the Ising interaction for spinsi^ 
The reason is that the constraint CI mentioned above 
and in Ref. [^, which states that a hole and a boson 
excitation are prohibited to occur simultaneously at the 
same site, cannot be violated here, because hopping t is 
strictly ID. This can be checked either by looking at Ht 
and verifying that the projection operators (1 — aja^) 
can be skipped without changing the physics, or by look- 
ing at the photoemission spectra in the limit of J — > 0. 
Whereas we did both of these checks, let us note here 
that for J = one obtains the incoherent spectrum with 
a bandwidth of Wine = 4-\/2< (not shown), which (unlike 
in the spin case) perfectly agrees with the RPA result 
Wine = 4V-Z — 2t from Ref. |4^, where z — 2 is the num- 
ber of possible forward going steps in the model. How- 
ever, still the three-site terms H^^^ and the orbiton terms 
Hj are not exact in Eq. (|5.8p and thus we checked the 
present results by comparing them with the numerical 
spectra obtained for the orbital Hubbard model (|5.ip — 
the results are presented in the next subsection. 

D. Comparison with numerical VGA results 

Since the problem of a hole added to the background 
with the AO order of t2g orbitals cannot be solved ex- 
actly using analytic methods and the SGBA had to be 
employed in the last section, we used also a numerical ap- 
proach. Actually, we compare the analytic results for the 
strong-coupling model (|5.2p presented in Sec. IV Gl with 
those obtained for the t2g Hubbard model (|5.ip using 
VGA. This enables us to compare not only the methods 
employed but also the two models which stand for the 
same physics in the strongly correlated regime. 

We first use the VGA to determine the staggered or- 
bital moment in the ground state of the t2g orbital Hub- 
bard model (|5.ip . 

mstagg = ^ E I - I ' (5.28) 

i 

with Q = (tt, tt) corresponding to the 2D AO order. 
We compare this result to the similar ones obtained for 
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FIG. 11: (Color online) Staggered magnetization mgtagg for 
increasing J/t as obtained for the t2g orbital Hubbard model, 
spin Hubbard model (called Hubbard on the figure), and for 
the Eg orbital Hubbard model, respectively. 



the spin Hubbard model and for the Cg orbital Hubbard 
model for a 2D plane (the strong-coupling model cor- 
rresponding to the latter situation was studied in Refs. 
[l^Ili.) As expected, "T'stagg increases with decreasing 
J (increasing U), see Fig. [TT] For J ^ (C/ oo), 
TOstagg 1 for both orbital models, corresponding to 
the perfect (classical) 2D Ising-like order, while quantum 
fluctuations reduce the moment mgtagg in the spin model. 
We remark that the treatment of the quantum fluctua- 
tions within the VGA is far from perfect (and limited 
by actual cluster size), so the staggered magnetization 
reaches Wstagg — 0.85 in the limit J — > (Fig. [TT|) . and 
does not reproduce the value of 0.606, well known from 
the spin- wave theory.— 

In all three models, the staggered moment (|5.28p ob- 
tained using the VGA decreases with decreasing U (in- 
creasing J), see Fig. [TTl because the kinetic energy can 
then generate more doubly occupied sites in the ground 
state. One finds that both the eg and t2g models give 
very similar results for the staggered moment, but dif- 
fer strongly from the SU(2) symmetric spin model, as 
has been shown before in three dimensions!^ However, 
we note that orbital order is slightly weaker for Cg or- 
bitals than for t2g. This may be easily explained by the 
fact that the Cg hopping is slightly smaller than t for the 
relevant orbital states l/-\/2(|z) ± while all other 
hopping processes are frustrated by the AO order (|5.28p . 
Gonsequently, correlations have a stronger impact on Cg 
electrons and induce a slightly enhanced mgtagg • Finally, 
we would like to emphasize that the AO is 2D in all 
three models, in spite of the fact that the kinetic energy 
is strongly anisotropic in the orbital models and actually 
has a ID nature in the t2g model, see below. 

Before we analyze the spectral functions, let us re- 
call that the VGA^ is appropriate for models with on- 
site interactions, as for instance the present Hubbard 
model (|5.1|) for t2g orbitals, but cannot be easily im- 
plemented for models where the interacting part con- 
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FIG. 12: Spectral function A{k,ijj) obtained with VGA for 
the 2D t2g Hubbard model (|5.1fl for: (a) a orbitals, and (b) b 
orbitals. Parameter: U = IQt. 



nects different sites, like in the t-J (or strong-coupling) 
model. For the present t2g model (|5.ip we use VGA 
with commonly used^ open boundary conditions, which 
leads to the spectral densities depicted in Fig. [121 The 
results resemble very much the SGBA results of Fig. [9] 
for the strong-coupling model (|5.2p . suggesting that not 
only both models are indeed equivalent in the strongly 
correlated regime, but also that the implemented SGBA 
method of Sec. IV Gl is of a very good quality. The dif- 
ferences between them, almost exclusively affecting high- 
energy features, are discussed below. 

On the one hand, we see that the high-energy part of 
the spectral density in Fig. [9] is composed of a number of 
peaks with a dispersion almost parallel to that of the QP 
state. In fact, the spectrum corresponds almost exactly 
to the ladder spectrum of the spin t-J model with Ising 
superexchange^^ but with a weak dispersion added to 
the peaks. The peaks at higher-energy are dispersive for 
the same reason as the QP state: After hopping a few 
times by NN hopping t — and creating string excitations, 
see Fig. [7] — the hole can exhibit coherent propagation 
via three-site terms, leading to the observed dispersion. 
On the other hand, the VGA spectrum (Fig. [T2) does not 
show these distinct peaks and the structure of ^(k, oj) is 
richer. However, the first moments calculated in separate 
intervals of w follow similar dispersions obtained for the 
first three peaks obtained in ^(k, w) within the SGBAJ^ 
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The above difference can be understood as following 
from the full Hilbert space used in the VGA calcula- 
tions which results in excitations of doubly occupied 
sites, weakening of the AO order even for relatively large 
U = lOi, see Fig. [iTl Therefore, the spectra of Fig. [12] 
have more incoherent features. In addition, the three-site 
terms which create two orbiton excitations (j5.6p . that 
were neglected in the SCBA, might also influence the 
high-energy part of the spectrum. The difference to the 
SCBA results might also be due to the fact that states 
with longer strings including several orbital excitations, 
which occur when the hole moves by a few steps via 
cannot be directly accommodated within the 10-site clus- 
ter solved here, and cannot be reproduced with sufficient 
accuracy. 

Apart from the differences in the high-energy part of 
the spectrum, we also observe differences in the spec- 
tral weight distribution (see also the detailed discussion 
below in Sec. IVE)) : In the VGA results (Fig. ^ the to- 
tal weight found in photoemission part (hole excitation) 
strongly depends on momentum k, while no such vari- 
ation can be seen in the SCBA results in Fig. [S] This 
difference does not originate from different approximate 
methods used, but stems from the different models: In 
Hubbard- like models, the number of electron states oc- 
cupied depends on the momentum ki^ In contrast, un- 
doped t- J-like models have exactly one electron per site, 
which enforces a different sum rule and eliminates the 
k-dependence from the photoemission part. 



E. Discussion of quasiparticle properties 




(b) 



0.8 



_ 0.6 



u 


= 20t 


— u 


= 12t 


u 


= 8t 


u 


= 4t 




= 2t 





In order to get a deeper understanding of the prob- 
lem mentioned in the last paragraph of Sec. IV D[ let us 
consider first the overall spectral weight distribution ob- 
tained in the VGA calculations. It is measured by the 
momentum-dependent electron occupation 

n„(k)^(4„Ck„), (5.29) 

obtained for the a-flavor in the Hubbard model, as for 
instance the t2g model (|5.ip . We recall that Eq. (|5.ip 
which leads in the limit [/ ^ t to the 2D t2g model 
(|5.2p is rather different from the one obtained for the spin 
Hubbard model with the SU(2) symmetry, see Fig. [T31 
One may easily identify the quasi- ID dependence only on 
kx in the h orbital momentum dependence nb(k) shown 
in Fig. [T^Ta). in contrast to the 2D variation of n„{k) in 
the spin case with isotropic hopping of Fig. I13f b). 

The insets in Fig. [13] show which states are occupied 
at ?7 = in the two models, and indicate the difference 
between the isotropic 2D hopping of the spin model and 
the ID kinetic energy of the orbital model. For instance, 

nff(k) = 0.5 along the (0, tt) (tt, 0) line for spins, while 

it shows full variation along this line in the orbital case. 
In both cases, we observe strong modifications of the elec- 
tron distribution with increasing U . For [/ = 0, the states 



FIG. 13: (Color online) Dependence of total weight found in 
the photoemission spectrum on momentum k for (a) a hole 
inserted into the h orbital of the t2g model (|5.ip and (b) for 
a hole with spin up in the SU(2) symmetric spin Hubbard 
model. All data were obtained by VGA. For illustration, we 
added a line at n{k) = 0.5, corresponding to the constant 
n{k) for f-J-like models. The insets show the first quadrant 
of the first BZ: The arrows indicate the path taken for the 
main panel, the shaded area gives momenta with nj, > 0.5 
and rij > 0.5, respectively. 



below the Fermi surface (k G S-p) are occupied and states 
above it (k ^ 5f) are empty. Consequently, ^^(k) is 
given by a step function with: nQ,(k) = 1 for k ^ S-p 
and nc((k) = for fc ^ iSp- The changes are particularly 
fast in the range of U ^ 8t; for U > 8t the momen- 
tum distribution function rict(k) (j5.29p smears out and 
one recognizes the strong-coupling regime. However, the 
difference between ^^(k = (0,0)) and na(k = (tt, tt)) is 
larger in the spin model, suggesting that the correlation 
effects are stronger in the orbital case. Indeed, this fol- 
lows from the ID character of the kinetic energy in the 
orbital model. In contrast, both strong-coupling models 
(for spin or orbital flavors) would give at half filling a 
constant UaCk) = 0.5 even for finite U < cx), although 



this result is strictly speaking correct only at [/ = oo. as 
shown in Fig. 1131 

After understanding the differences between the QP 
properties found in the VGA and the SCBA, we concen- 
trate solely on the QP properties calculated using the 
latter method. Hence, following Ref. H, we analyze the 
characteristic features of the QP states in the 2D t2g 
model, such as the bandwidth W and the QP spectral 
weight aqp . The energy of incoherent excitations (string 
states) is to some extent characterized by the separa- 
tion between the QP state and the next (second) spectral 
feature at higher energy - it is called here a pseudogap 
A. All these quantities increase with increasing superex- 
changc energy J {t — J/i), sec Fig. [T31 One finds that: 
(i) the bandwidth Wi of the first QP peak, see Fig. [TiTa). 
is proportional to for small J {J < 0.7) and to J in the 
regime of large J ( J > 0.7) — the bandwidth renormal- 
ization is here distinct from the one found either in the 
spm SU(2) (see Ref. H) or in the orbital eg models,— (ii) 
the bandwidth W2 of the second largest dispersive peak 
[Fig. [HTa)] is smaller than that for the first peak and 
tends to saturate at W2 ^ 0.25t value for larger J > t 
(not shown), (iii) the spectral weight aqp of the QP peak, 
shown in Fig. I14r b). grows with J. and (iv) the pseudo- 
gap A shown in Fig. [TUc) grows generally like J^^^, 
while for higher J values some deviation from this law 
is observed for the k = (0,0) point. Most (but not all) 
of these results are qualitatively different from the ones 
obtained for the QP states, and their momentum depen- 
dence, in the SU(2) Heisenberg antifcrromagnet. Let us 
now discuss the above mentioned QP properties in more 
detail. 

Firstly, the QP bandwidth arising from the superex- 
change three-site terms is renormalized as it is much 
smaller than the respective free value, W ^ 2J. Even 
at J = t, the QP bandwidth is only W ~ J/2, i.e., is 
here reduced by a factor of 4. This is not surprising in 
the view of incoherent processes which " dress" the prop- 
agating hole and increase its effective mass. Indeed, the 
collapse of the QP bandwidth in the regime of J ^ 
may be understood as following from numerous incoher- 
ent string excitations which are easy in this regime as 
they do not cost much energy. A similar but consid- 
erably weaker reduction of the ID dispersion by string 
excitations was seen before in the centipede model, see 
Fig. ini^a). However, in that case the renormalization was 
almost linear as the length of the string excitations was 
limited to a single step (within one of the three-atom 
units along the chain), and could not further increase 
with decreasing J. In addition, the dispersion of the sec- 
ond peak is weaker than that of the QP. Interestingly, the 
bandwidth corresponding to the dispersion of the second 
peak in the centipede model is not only weaker than that 
of the QP itself, but is also renormalized in a similar 
way to that found for the full 2D t2g model. Altogether, 
this suggests that the bandwidth renormalization of the 
coherent hole propagation in the 2D t2g strong-coupling 
model follows from the creation of string states during the 




J/t 

FIG. 14: (Color online) Quasiparticle properties obtained for 
the 2D t2g model within the SCBA for increasing superex- 
change J (with r = J/4): (a) the bandwidth of the QP Wi 
(solid line) and the second dispersive feature W2 (dotted line), 
(b) the spectral weight aqp, and (c) the distance between 
the first two peaks in the spectra (pseudogap) A. The solid 
(dashed) lines in (b) and (c) give the results for k = (0, 0) 
[k = (7r/2, 7r/2)], respectively. The light solid line in (c) indi- 
cates the t{J/t)^^^ law (see text). 



ID hole propagation via three-site terms. Such processes 
are absent in the ID and 2D FK model, and therefore 
the hole moves there freely by three-site hopping terms 
and the bandwidth is unrenormalized. 

Secondly, in contrast to the spin t-J model with 
Ising superexchange interactions^ where the QP spectral 
weight is independent of k, it varies here with the com- 
ponent {kx or ky) of the momentum k, cf. Fig. [T3]as well 
as Fig. [51 and [T^ Similar to the spin t-J model, the QP 
spectral weight is larger for the k values with the lowest 
QP energies than for the ones close to the maximum in 
QP dispersion, for instance aQp(7r/2, 7r/2) > aQp(0,0), 
see Fig. [HT b). Altogether, the k-dependence here is 
however much weaker than in the spin case.— The in- 
crease of aQp(7r/2, 7r/2) with J/t resembles the increase 
of the spectral weight for the low-energy peak at /c = 7r/2 
in the centipede model, see Fig. [IJb). 

Finally, we address the issue of the pseudogap which 
separates the QP state from the first incoherent excita- 
tion. It scales almost as i(J/t)^/^, see Fig. [HTch in 
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agreement with the result for the Ising spin model This 
demonstrates that in spite of the observed k-dependence 
of the QP properties and the pseudogap itself, the pseu- 
dogap originates from string excitations similar to those 
generated by the hole moving in the spin background 
with AF order. Note also that the spectrum of the 2D 
model with dense distribution of incoherent maxima in 
the range of J is qualitatively different from the ID 
centipede model, shown in Fig. [5l^b). 



VI. PHOTOEMISSION SPECTRA OF 
VANADATES AND FLUORIDES 

In this Section we discuss the possible implications of 
the results obtained for the orbital model of Sec. |V] 
on future experiments, and make predictions concerning 
the photoemission spectra of strongly correlated fluorides 
and vanadates. As before, we discuss the strongly cor- 
related regime with U ^ t. The first important fea- 
ture to consider is the interplay of the three-site hop- 
ping with the longer-range {^2,^3} hopping to second and 
third neighbors which contributes to the electronic struc- 
ture and may always be expected in any realistic system 
(for instance, due to hybridization with oxygen orbitals). 
These hopping elements were neglected in both the Hub- 
bard model (|5.1|) and in the strong-coupling model (|5.2p . 
but they could significantly influence the spectral weight 
distribution. We will see, however, that although fea- 
tures induced by longer-range hopping arc small as long 

1^2(3)1 < t, they can be clearly distinguished from the 
effects of three-site hopping. 

The same requirements for orbital symmetry that are 
necessary to obtain NN hopping, as discussed in this 
work, also strongly restrict the range of allowed longer- 
range hopping terms. It is important to recall that the 
d-d hopping elements involve intermediate oxygen or- 
bitals. For next nearest neighbor (NNN) hopping, the 
orbital phases of the involved oxygen 2p^ orbitals make 
all terms vanish that conserve orbital flavor^i^ and only 
orbital-flipping terms, 

i 

given by hopping element t2, are finite. With realistic 
parameters, we arrived at the estimation of \t2\ ~ 20 
meV, i.e., \t2\ ~ J/3. Similar to the orbital flipping three- 
site term (|5.6p . such a hopping process disturbs the AO 
order stabilized by the supcrcxchangc and induces string 
excitations. For this reason, its impact is largely confined 
to the high-energy part of the spectrum and is rather 
small for the low-energy QP state. This can be seen in 
Fig. [THl where we show the spectral density for t2 = O.lSt 
and J = OAt: While the higher energy part is somewhat 
affected by finite <2, the intensity and dispersion of the 
low-energy QP is almost the same as obtained for t2 = 0, 
see Fig. \Mh). 
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FIG. 15: Spectral density A(,(k,(jj) obtained within the VGA 
method for a hole inserted into b orbitals of the t2g model 
(|5.ip . supplemented by finite NNN hopping (|6.ip . Parame- 
ters: U = Wt, and t2 = 0.15t. 

The QP dispersion could also be influenced by the 
third-neighbor hopping terms t^, where the orbital sym- 
metry leads to the same anisotropy as for NN hopping: 
a orbitals allow only hopping along the a axis, and b or- 
bitals only along the b one: 

{imj}\\a {imj}\\b 

Here the unit consisting of three sites {imj}, shown in 
Fig. EJa), is parallel to one of the cubic axes in the (a, b) 
plane. In contrast to t2 terms, these terms do not induce 
any string excitations but contribute only to the QP state 
itself, so they mix with the three-site effective hopping 
T. To illustrate this effect, we have chosen = ±J/4 
for the spectra shown in Fig. [111 Note that the value 
of 1^3 1 is here larger than expected in transition metal 
oxides, where it is in general smaller than the three-site 
hopping term r = J/4. The spectral density ^(k, w) 
contains now the combined effects of the three-site terms 
cx T and third- neighbor hopping cx ts, and one finds that 
^3, depending on its sign, can cither amplify or weaken 
the QP dispersion which stems from the effective three- 
site hopping, see Fig. [16] 

From the above example we have seen that the longer- 
range hopping violates the particle-hole symmetry of the 
spectral functions. The spectra obtained for the original 
orbital Hubbard model (|5.ip with NN hopping t obey the 
particle-hole symmetry. The three-site supcrexchange 
terms arise from this model, and therefore these terms 
also have to follow the particlc-holc symmetry. This is 
in marked contrast to the t2 terms that do not respect 
it^ or to terms, see Fig. [1^1 As a result, the spec- 
tra exhibit a striking particle-hole asymmetry — reduced 
dispersion in the particle (inverse photoemission) sector 
corresponds to enhanced dispersion in the hole (photoe- 
mission) sector and vice versa. 

We will show now that the above asymmetry follows 
indeed from the difference between the NN and NNN 
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a axis then becomes 



FIG. 16: Photoemission [{uj — /i) < 0] and inverse photoe- 
mission [(lo — /i) > 0] part of the spectral density Ab{k, lo) for 
a hole inserted into b orbitals, obtained within VGA for the 
t2g model (|5.ip with an additional longer-range third-neighbor 



hopping ta (|6.2|l . The value was selected to suppress disper- 
sion arising from the three-site effective hopping (|5.5|l in; (a) 
the hole (photoemission) sector with to, = O.lt = J/4, and (b) 
in the inverse photoemission sector with t-j, — —O.lt = —J/4. 
Parameter: U = lOt. 



hopping under particle-hole transformation. While this 
is transparent for the Hubbard model acting in the full 
Hilbert space, it is somewhat subtle for the t- J-like mod- 
els. Thereby we focus on the hopping which influences 
directly the QP dispersion. The operator for NN hop- 
ping can be transformed from {c^., cj} electron operators 
to {hl,hr} hole operators, and one arrives at an identi- 
cal form for the kinetic energy, as long as a phase shift 
between the two sublattices is introduced: 



K = i-i] 



= (-l)(''-+'^« 



(6.3) 



^(c|;Cr+a + 4+aCr) 
r 

r 

r 

^ ( /^t/v+a + hl^^K ) , (6.4) 



and analogously along the b axis. The minus sign for 
one of the sublattices corresponds to a momentum shift 
by q = (tt, tt), as can be easily verified in the Fourier 
transform. 

1 



hi = — Vc*''(-i)('^=-+''«)^ 



N ^ 

r 

N ^ 



(k+q)r _ 



(6.5) 



where r = (r^^ry) is the lattice site. Hopping along the 



The on-site density-density interaction is not affected by 
the particle-hole transformation, apart from a shift in the 
chemical potential. 

Since the three-site hopping emerges from the 
Hubbard-like model with NN hopping, it respects 
particle-hole symmetry. Hence it obeys the same rules 
concerning particle- hole transformation, i.e., momentum 
(0,0) for electrons is mapped to (7r,7r) for holes. For 
the third-neighbor hopping (|6.2[) . however, the above 
transformation does not longer work, because both the 
creation and the annihilation operator act on the same 
sublattice. Instead the transformation vector would have 
to be q' — {tt/2, 7r/2). Consequently, the combined effect 
of explicit NNN hopping and three-site terms stemming 
from NN processes turns out to be strongly particle-hole 
asymmetric. For example, negative ^3 gives a band in 
the electron sector with the largest distance from the 
Fermi energy at momenta (0,0) and (tt, tt), and the val- 
ues nearest to it at (7r/2,7r/2), and the same is true for 
the three-site hopping. Consequently, the two dispersions 
add together and lead to increased total dispersion, see 
the photoemission part in Fig. |16(b)[ On the contrary, in 
inverse photoemission the direct NNN hopping ^3 gives 
a maximal distance at (7r/2,7r/2), while maximal energy 
is still found at (0, 0) and (tt, tt) for the three-site terms. 
Therefore, now and three-site hopping t compete with 
each other, and the dispersion is weaker. For a particu- 
lar choice of the model parameters they can even cancel 
each other, as shown in the inverse photoemission part 
in Fig. 16(b) Positive ^3 leads to the opposite result, see 
Fig. |16(a)| Thus, even large and unphysical values of ts 
not only do not destroy the qualitative spectra predicted 
in the previous sections but result in asymmetry between 
the photoemission and inverse photoemission part of the 
spectra, so their contribution can easily be resolved. 
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The symmetry arguments leading to Eq. (|6.ip and Eq. 
(|6.2p remain valid also for systems with specific Cg or- 
bital degeneracy, as observed in certain fluorides with 2D 
AO order which involves alternating — and — 
orbitalsi^ In fact, the orbital model given by Eq. (|5.2[) 
describes also this case, as we show by a detailed deriva- 
tion in the Appendix. Hence, we conclude that the pho- 
tocmission and inverse photoemission spectra for the pla- 
nar vanadium oxide Sr2V04 and for the planar K2CUF4 
or Cs2AgF4 fluorides should be qualitatively similar to 
the spectral functions shown in Figs. [SlorfT^ 



VII. SUMMARY AND CONCLUSIONS 

In this paper we analyzed only the orbital sector of 
the superexchangc, which decides about the hole dynam- 
ics when spins are polarized in the FM ground state. 
We discussed all possible situations (see below) where 
the orbital symmetry leads to the purely Ising superex- 
changc in one and two dimensions. Exceptions from this 
rule are numerous systems with Cg orbital degrees of 
freedomjiiSii^ or FM chains with two active orbitalsf^i^ 
but we also provided examples of Cg systems with Ising 
superexchangc. 

The ID Hubbard-likc model with two orbital flavors, 
but only one of them participating in NN hopping, served 
to explain the general principles and consequences of the 
Ising-like superexchangc. Besides, this model stands for 
several physically relevant situations, including electrons 
moving within cither Cg or t2g orbitals in one dimension, 
and the ID FK model. We have shown that, particu- 
larly in all these cases, the relevant strong-coupling model 
has to include the three-site effective hopping. When 
both interorbital hopping and orbital-flip processes in the 
superexchangc are absent, the three-site hopping term 
which arises from superexchangc is crucial and is the only 
source of coherent hole propagation. 

We have shown that the 2D FK model with one im- 
mobile (/) and one mobile (d) orbital has many common 
features with the ID model. In both cases one finds only 
one dispersive mode for a hole inserted into the mobile 
orbital, and two non-dispersive modes for a hole doped in 
the immobile orbital. This latter hole excitation creates 
a trapped polaron, with the hole confined within a clus- 
ter consisting of a central site and its nearest neighbors 
(i.e., three sites in the ID model, and five sites in the 2D 
model). While the hole can in principle escape from the 
polaron via three-site hopping process, we have shown 
that such processes have only very low spectral weight 
in the realistic regime of parameters, and thus the hole 
remains de facto trapped inside the polaron. In contrast 
to this almost perfectly localized hole, a hole in the mo- 
bile orbital propagates freely, and its dispersion which 
originates from the three-site hopping is unrenormalized. 
Therefore, the two inequivalent orbital flavors behave in 
the FK models in a radically different way, and decouple 
from each other (interacting only by the on-site Coulomb 



interaction U , which stabilizes the AO ground state). 

The model relevant for the 2D orbital physics in tran- 
sition metal oxides leads, however, to qualitatively dif- 
ferent results. In the 2D t2g orbital model, which is also 
applicable to the AO state formed by Cg orbitals in fluo- 
rides (see the Appendix), electrons do not separate into 
those conflned to either sublattice (occupied by orbitals 
of particular symmetry in the ground state with AO or- 
der), but may dclocalize over the lattice and thereby 
undergo incoherent scattering on the orbital excitations, 
which strongly renormalizes and reduces the dispersion 
of the QP states. These QP states arise at half filling in 
the regime of large Coulomb interaction f/, when a hole 
(electon) is added to the ground state with AO order. 
While electron hopping is of the purely ID character, it 
selects by symmetry possible three-site processes, which 
are responsible for the QP dispersion. Therefore, the dis- 
persion of the QP state depends on the considered orbital 
and is again ID, with a hole propagating coherently along 
the two crystal axes for the two orbitals. 

We emphasize that the mechanism of coherent hole 
propagation which occurs in the 2D t2g orbital model 
is completely different from the one known in the spin 
case. In orbital systems (with conserved orbital flavors) it 
originates entirely from the three-site hopping processes, 
similar to hole propagation in the ID or 2D FK model. 
But unlike in the latter models, in the 2D t2g case the 
QP bandwidth is strongly reduced from the value given 
by the amplitude of bare three-site hopping. We have ex- 
plained this renormalization as following from incoherent 
string excitations which dress the coherent propagation 
and do not contribute additional momentum dependence. 
As a special case, we have discussed the subtle interplay 
between the coherent hole propagation and string exci- 
tations in the ID centipede model, where polaronic hole 
confinement competes with coherent propagation along 
the chain, and which to some extent resembles the real- 
istic 2D t2g case. 

We discussed the impact of realistic longer-range hop- 
ping terms (as expected in real materials such as vana- 
dates or fluorides), and found that the second neighbor 
terms are frustrated in the ground state with AO order 
- these processes would flip the orbital flavor, and are 
therefore suppressed at low energy, not affecting the QP 
dispersion. In contrast, third-neighbor hopping processes 
conserve orbital flavor and lead to a pronounced particle- 
hole asymmetry in the spectral weight distribution. In 
both cases, the ID character of hole propagation which 
follows from the symmetry of involved orbitals, survives 
and determines the character of the spectral density at 
low energy. 

In summary, we have demonstrated that orbital mod- 
els with Ising superexchangc describe a broad class of 
interesting phenomena. Spectral features resulting from 
such models exhibit weak momentum dependence and are 
fundamentally different from those known from the the 
spin case with the SU(2)-symmctric superexchangc. The 
predictions of the theory presented in this paper provide 
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an experimental challenge for the transition metal oxides 
with orbital degrees of freedom, where similar features 
could possibly be observed in FM planes with AO order. 
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APPENDIX: STRONG-COUPLING MODEL FOR 
FLUORIDES 

Here we show that the model developed in Sec. FVlmav 
also be applied to certain fluorides with FM planes and 
AO order. In contrast to the orbitally degenerate sys- 
tems, in the systems with Cg orbital degeneracy the lat- 
tice distortions in the cubic phases are usually quite large. 
In particular, the static distortions may counteract to 
some extent the AO order favored by the superexchange 
interactions, as e.g. in undoped manganites RMnO^r- or 
fluorides Cs2AgF4i^ However, the crystal fleld does not 
suppress the orbital order present in these systems but 
instead it only modifles the occupied orbitals which form 
the AO state. They have to be optimized in a micro- 
scopic model by choosing particular linear combinations 
of the Eg orbitals, which form the AO order, in order to 
fit best to the superposition of the superexchange and the 
Jahn- Teller terms generated by ligand fields 4^ In certain 
situations this "modification" could be quite substantial 
and could even lead to such a selection of such Cg orbitals 
that the resulting state is modified to a ferro-type orbital 
order 

At finite crystal field splitting cx E^, it is convenient 
to describe the changes in the occupied orbital states 
by making two complementary transformations at both 
sublatticeSfi rotating the orbitals by an angle 9 = j — (p 
on sublattice A, and by an angle 6 = j + (j) on sublat- 
tice B, so that the relative angle between the occupied 
orbitals is ^ — 20 and decreases with increasing 0, i.e., 
with increasing Ez, 

f\^l)A_f cos(f-0) sin(f - 0) W |z) A 
l^-sin(f-0) cos(f - 0) ; 1^ |x) J ' 

(A.l) 

{\fi)j\_f cos(f + 0) sin(| + 0) \ / |z), \ 
V J [- sin(f + 0) cos(f + cf>)j[ \x), J ' 

(A.2) 

where the "old" orthogonal (basis) orbitals are defined 
as |x)i = 75 1^;^ - y'^)t and |z)j = 7g|3z^ - r^)i for every 
sublattice site i. Due to the above transformation the AO 
order is formed now by \fi)i and occupied orbitals at 



sublattices, i G A and j G B, respectively. Let us stress 
that although the transformation defined by Eqs. (|A.H - 
IA.2P is orthogonal this does not mean that orbitals on 
different sublattices, such as e.g. the occupied orbitals 
and \v)j, are orthogonal for any arbitrary angle 0. 
For the 2D FM systems with active eg orbitals which 
are considered here, the relation between the crystal field 
Ez and the optimal orbital configuration defined by the 
angle [see Eqs. (12) and (13) of Ref. Q is given by: 

i;^=4Jsin20, (A.3) 

where J is the superexchange constant. In the case of 
fluorides such as Cs2AgF4 (Ref. or K2CUF4 (Ref. 
[28[ ) discussed here, the filling is one eg electron per site 
and the crystal field would select the angle = 7r/12 (for 
the reason of looking at this angle see below) since the 
convenient basis adapted to the actual AO order looks as 
follows: 




(A.4) 

where the occupied (empty) orbitals for this type of AO 
order are denoted as |a;) (|z)) on both sublattices. 

The reason why these particular pairs of basis orbitals 
(|A.4[) are interesting here is that this is the only choice 
of occupied Cg flavors which forms a two-sublattice AO 
order with the interorbital hopping between occupied or- 
bitals vanishing by symmetry, and where the interactions 
described by pseudospin operators do not allow for any 
quantum fluctuations. This resembles the t2g case dis- 
cussed in this paper. There is, however, one subtle differ- 
ence: two occupied {|a;)i, \x)j} orbitals on sublattices A 
and B are not orthogonal and do not form the global ba- 
sis in the Cg orbital space. The choice made in Eq. (|A.4p 
means that we consider two different pairs of orbitals for 
both sublattices, and the interorbital hopping between 
the unoccupied orbitals is also rather small but remains 
finitei^ Hence, the respective strong-coupling Hamilto- 
nian is richer than the one for the t2g case and we need 
to check under which conditions it can be reduced to a 
similar polaron Hamiltonian as Eq. (jS.Sp . 

The eg orbital t-J Hamiltonian for the FM planes with- 
out the three-site terms but including the crystal field was 
given e.g. in Ref. [60I . Here we rewrite the kinetic term in 
a slightly different form (there it was written already us- 
ing the polaron representation) and substitute — 7r/12 
to obtain: 

He, =Ht + Hj + Hz , (A.5) 
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where 



i 

V3 



i eA 



V3 



fe>ll'i 



+ 7:J 



2 ^ 



(A.7) 



ieB 

Here TIf = i(fij^ - riij;) for i e A, = i(njj 



(A., 



j & B, and T!f = for every site i, see Ref. 0- 

As before, a tilde above a fermion operator indicates that 
the Hilbert space is restricted to unoccupied and singly 
occupied sites, e.g. xl = a;J(l — riiz). The last term 
Hz represents the above mentioned crystal field with the 
strength of the interaction written according to Eq. (|A.3P 
with (j) = 7r/12. 

However, we are not aware of any work where the 
three-site terms complementing such a t-J model have 
been derived. We use second again order perturbation 
theoryi^i^^ applied to the Hubbard model for spinless 
Cg electrons in a FM plane^^ with the basis rotated by 
(I) = 7r/12, following Eqs. (12) and (13) of Ref. 0. This 
leads to the following three-site terms for the eg strong- 
coupling model (with cj) = 7r/12): 



Hr = + H° + Hf, 



for where 



(A.9) 



(A.IO) 



;nixZ,^i^ + 3z^. ,hizZ,^i- \/?,z^, - V3z|_g5jzjZ._j^g + h.c. 



(A.ll) 



ieA 



— r [z^-u-^iir^-Li; — 3zh~z^-XiX_,i + -s/Sz^.-r 



+ ^a^-iaj^jzrb 3z|_|_gz|x,a;- |j + V3z|_|_gri„x^^ \/izl^^z\xiZ^ 4- h 



(A.12) 



Here we underlined (doubly underlined) terms which do 
not require orbital excitations (require orbital excita- 
tions), respectively, i.e. 

H^^^^H^, Hi^'>=H^. (A.13) 

Next, we perform the same standard transformation to 



obtain the polaron Hamiltonian from the strong-coupling 
modcl^ for the lightly doped ordered states as done in 
Sec. IV C[ i.e. we introduce boson operators ai (orbitons) 
and fermion operators hi (holons) which are related to 
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the Xi and Zi operators in the following way: 

x\ = h,{l - a\ai), zl = h,a\. (A.14) 

Please note, however, that here we did not have to per- 
form rotation of the pseudospins since we defined distinct 
electron operators for the occupied and empty orbitals, 
cf. Eq. jXH. 

Again we implement a linear orbital-wave 
approximation^ (we keep only linear terms in or- 
biton operators) and we skip (\ — h\hi) operators (which 
in fact is not an approximation if there is only one hole 
in the entire plane). This means that e.g. the three-site 
terms are reduced only to the terms which were either 
underlined or doubly underlined in Eqs. (jA.10[) - (|A.12[) . 

i.e. to either H'f'' or Here, however, we have to 

use yet another approximation which was unnecessary 
for the t2g model: we skip terms h'^^ which were absent 
in Sec. IV Bl This approximation is allowed since these 
terms contribute to the vertex as cx r and not as oc 
resulting typically in much reduced energy scale for the 
new vertex contributions. Furthermore, we showed in 
Sec. Ell that such terms [cf. Eq. and Fig. [15] do 

not change the energy of the QP and merely modify 
the incoherent spectrum. Eventually, we arrive at the 
polaron Hamiltonian for the holes doped into the eg 
orbitals of the fluorides, with the hopping terms: 

Ht = ^^-7^ XI {cos(fc^ - qx)hiAK-c^M(^c^A 

k,q 

+ cos(A:y-gy)/ij,B/ik-q,AaqB + h.c.| , (A.15) 

^r"^ = ^^Il{cOs(2fcj,)4^/lk^+COs(2fc,)/lj^5/lkB} , 

k 

(A.16) 



and the remaining terms resulting in the energy renor- 
malization 



H.J + H:, = 4-^X1 ("L"kA + aL^kij) • (A.17) 

k 



Therefore, the Hamiltonian given by Eqs. (jA.15|) - (jA.17p 
reduces to the polaron Hamiltonian (|5.8p after substi- 
tuting -\/3i/2 t, and consequently 3 J/4 J and 
3t/4 T. This substitution stems from the different 
definitions of the hopping t in the Cg and in the t2g sys- 
tems — in the former case it is the {dda) hopping between 
the — orbitals along the c direction, whereas in the 
latter case it is the hopping clement between a pair of 
active t2g orbitals, e.g. yz orbitals in the (a, b) plane. 



In summary, we have shown that the Hamiltonian 
given by Eqs. (|A.15|) - (|A.17p provides the framework to 
analyze the behavior of certain lightly doped eg systems, 
with FM planes and AO order which suppresses the in- 
terorbital hopping between occupied orbitals. Its equiv- 
alence to the polaron model (|5.8p demonstrates that the 
results obtained and discussed in Sec. |V] should also ap- 
ply to the case of a hole doped into the fluoride plane 
with the AO ordered eg orbitals. 
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